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ABSTRACT 

We introduce the Mass-Excitation (MEx) diagnostic to identify active galactic nuclei (AGNs) in 
galaxies at intermediate redshift. In the absence of near- infrared spectroscopy, necessary to use 
traditional nebular line diagrams at z > 0.4, we demonstrate that combining [0 ill] A5007/H,S and 
stellar mass successfully distinguishes between star formation and AGN emission. The MEx classifi- 
cation scheme relies on a novel probabilistic approach splitting galaxies into sub-categories with more 
confidence than alternative high-z diagnostic diagrams. It recognizes that galaxies near empirical 
boundaries on traditional diagrams have an uncertain classification and thus a non-zero probability of 
belonging to more than one category. An outcome of this work is a system of statistical weights that 
can be used to compute global properties of galaxy samples. We apply the MEx diagram to 2,812 
galaxies at 0.3 < z < 1 in the Great Observatories Origins Deep Survey North and Extended Groth 
Strip fields, and compare it to an independent X-ray classification scheme. We identify Compton-thick 
AGN candidates with large X-ray absorption, which we infer from the luminosity ratio between hard 
X-ray emission and [O in] A5007, a nearly isotropic tracer of AGN. X-ray stacking of sources that 
were not detected individually supports the validity of the MEx diagram and yields a very flat spec- 
tral slope for the Compton-thick candidates (F « 0.4, unambiguously indicating absorbed AGN). We 
present evidence that composite galaxies, which are difficult to identify with alternative high-redshift 
diagrams, host the majority of the highly-absorbed AGN. Our findings suggest that the interstellar 
medium of the host galaxy provides significant absorption in addition to the torus invoked in AGN 
unified models. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: fundamental parameters and 
ISM — galaxies: high-redshift — X-rays: galaxies 



1. INTRODUCTION 

Most if not all galaxies contai n a supermassive blac k 
hole (SMBH) in their center (jRichstone et all Il998h . 
which follow fundamen t al relations to their host g alaxies 
(|Magorrian et alj|1998t [Ferrarese fc Merrittl 120001) . Fur- 
thermore, the accretion history of supermassive black 
holes follow a similar, scal ed-down, trend as the cos- 
mic star formation history (Ba rger et al.l l200lT ) . Taken 
together, these observations strongly suggest that black 
hole accretion and star formation may be linked phe- 
nomena. In particular, active galactic nuclei (AGNs) 
have been invoked in galaxy evolution models as a 
means to control t he rate of star formation in galaxies 
via fe edback (e.g., ICroton et al.l 120061 : INaravanan et al.l 
2008). While some observational evidence is provided 
by, e.g., high- velocity outflows in post-starburst galaxies 
(|Tremonti et al.|[2007l ). the full picture on the interplay 
between AGN and their host galaxies remains unclear. 

A complete understanding of galaxy evolution re- 



quires the study of both galaxy stellar content and nu- 
clear activity. However, differentiating the powering 
source - star formation versus AGN - poses an ap- 
preciable challenge. While there are several indepen- 
dent tracers of AGN, all suffer from limitations. Be- 
cause the caveats of one tracer can be overcome by the 
strength of another, it has become clear that multi- 
wavelength methods are required to assess the ubiquity 
of AGNs. For example, X-ray emission has been used 
extensi vely to uncover and study large populations of 
AGN s (|Bade et al.lll995l iBovle et al.lfl993t iBrandt et all 
l2001h . However, even the most sensitive X-ray sur- 
veys still miss heavily absorbed systems. These ab- 
sorbed AGNs were inferred by the unresolved portion 
of the cosmic X-ray background, which shows a flat 
X-ray sp ectral slope highly suggestive of X-ray ab- 
sorption dComastri et all 119951: iMushotzkv et aU boOfi 
I Alexander et al.ll2003l : lBauer et al.l2004HTreister fc Urrvl 
I2005D . 
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This scenario is supported by the existence of 
Compton-thick systems in the nearby universe that 
would not be detected in the most sensitive X-ray sur- 
veys if they were at higher redshift, e.g., NGC 1068. This 
galaxy was also shown to have X-ray absorbe rs on a scale 
of the order of ~1 pc (jGuainazzi et al.H2000() . 

Even if X-ray photons are absorbed or scattered by 
material on a compact scale such as the torus compo - 
nent described in the unified model (jAntonuccil Il993). 
emission originating from larger scales may reach the 
observer regardless of the line of sight. Such isotropic 
tracers include emission lines from the narrow line re- 
gions, which are exterior to the torus. C ommonly used 
lines i nclude [O III] A5007 in t he optical (Baldwi n et al.l 
[1981 IBusko fc Steinerl fl988h and [OlvlA25.9um or 
[Ne vl A14/im in th e infra r ed regime (ISturm et al.l 1 20021 : 
Armus et al l 12004 [20071: IDiamond-Stanic et all 12009; 
Rigbv et all l2009t iLaMassa et al.l 120101) . Another such 
tracer is mid-infrared continuum emission originating 
from hot dust heated by the AGN (e. g., lLacv et all 
[200llStern et all2005HDonlev et all2007t) . This method 
works especially well for intrinsically luminous AGN that 
are deeply enshrouded, but does not select the less- 
luminous systems which may still be absorbed. Mid- 
infrared aromatic features can also be used to diag- 
nose the powering source in galax ies (|Genzel et allll998l : 
iLutz et alJU99l [Pope et all 12001 ). 

On the one hand, some of the distinctions between the 
classes of AGNs and the selection methods listed above 
are thought to arise from their orientation with respect 
to the observer's line of sight. On the other hand, some 
AGN classes seem to comprise physically distinct phases 
of activity (i.e., low versus high accretion rate, beginning 
or end of an active phase). The latter are more interest- 
ing from the point-of-view of galaxy evolution as they 
may provide a handle on the importance of AGN phases, 
their duty cycle, and the interplay between AGNs and 
their host galaxies. 

One example of physically-distinct AGN phases are 
Seyferts and low-ionization nuclear emission line regions 
(LINERs). These two categories emerged from optical 
spectroscopy studies where emission lines with different 
excitation properties were used as probes of the radia- 
tion e xciting the interstellar gas (Scvfcrt 1943: iHeckmanl 
1980). In addition to exhibiting specific spectral signa- 
tures, Seyfert and LINER nuclei were found to reside in 
distinct host galaxies compared to st ar-forming galax- 
ies an d also relative to one another. iKauffmann et all 
( 2003) showed that galaxies hosting an AGN tend to have 
a larger stellar mass compared to star-forming galax- 
ies that lack optical AGN signatures. Galaxies with a 
Seyfert nucleus (or Seyferts for short) often have a young 
or intermediate- age component in their stellar population 
whereas galaxies with LINER emission have, on aver- 
age, an older stellar popu lation as well as a la rger stellar 
mass than Seyferts (e.g. iKewlev et all [20061 Putting 
this evidence together with the observations that AGNs 
seem to follow a decreasing sequence in accr etion rate 
from Seyferts to LINERs to composite galaxies (|Hdl2008l ) 
suggests an evolutionary picture where LINERs may be 
older, dying, AGNs relative to Seyferts. 

Whether this emerging picture is supported by higher- 
redshift observations is unclear. At larger distances, 
it is generally more difficult to get a complete cen- 



sus of galaxies with AGNs, let alone to classify 
them in different AGN categories/phases. Gather- 
ing complete samples of galaxies for which we know 
the powering source, and whether the central black 
hole is actively accreting, is especially challenging at 
z > 0.4. Beyond that redshift, optical emission 
lines needed for AGN classification such as Ha and 
[N ii] A6584 are shifted into the near-infrared, preventing 
the application of well-ca l ibrated, traditional diagnos- 
tics ([Baldwin et al 1 119811: iVeilleux fc Osterbrockl 119871: 
Kewlev et al.ll2001t IKauffmann et alJl2003c IKewlev et all 
20061 : IStasihska et al.ll2006l) ~ 

In this paper, we present the Mass-Excitation (MEx) 
diagnostic diagram, based in part on optical nebular lines 
that can readily be observed out to z ~ 1. Following a 
similar method as Wei ner et all (|2007l ). who replaced the 
[N ii] A6584/Hq line flux ratio used in the BPT diagram 
(|Baldwiri~e t al. 198l|) with absolute -ff-band magnitude, 
we adopt stellar mass as a substitute for [N n] A6584/Ha. 
We will show that a better census of AGN can be ob- 
tained by finding both intrinsically weak AGNs as well 
as absorbed systems that are undetected in X-ray ob- 
servations. Our classification scheme relies on a novel 
probabilistic approach and allows us to split the galaxies 
into the following categories: purely star-forming galax- 
ies, Seyfert 2s, LINERs or composite systems (i.e., with 
both star- formation and AGN) . Galaxies near the empir- 
ical boundaries on traditional diagnostic diagrams have 
a less certain classification and are thus assigned a non- 
zero probability of belonging to more than one category. 
As a result, the MEx diagnostic also outputs statistical 
weights that can be utilized to compute global properties 
(e.g., stellar mass, metallicity, etc.) in statistical samples 
of galaxies belonging to any of the categories listed above. 

Using a sample of low-redshift galaxies described in 
^21 we calibrate our diagnostic in §3.11 We briefly an- 
alyze the occurrence of low-ionization nuclear emission- 
line regions (LINERs) in i)3. 21 before introducing a novel 
approach to galaxy spectral classification based on the 
probability of each spectral class (e.g., star- forming or 
AGN, §373"]) . We describe our 0.3 < z < 1 galaxy 
sample in §31 and the application of our new Mass- 
Excitation (MEx) diagram in §5] We find an excellent 
agreement between the MEx diagram and the X-ray clas- 
sification f £|5.1j) . We examine the different optical and 
X-ray classes more closely with an X-ray stacking anal- 
ysis ( §5.2j) . Combining hard X-ray emission and opti- 
cal emission lines allow us to probe the X-ray absorp- 
tion leading to the discovery of Compton-thick AGNs 
among our intermediate-redshift sample §21 We com- 
pare the MEx diagram with alternative AGN diagnostic 
diagrams in §7.11 and we describe how the new method 
contributes to achieving a more complete census of AGNs 
( §7.2|) . We discuss the fraction of AGNs that suffer from 
X-ray absorption (fc |7.4[) and possible evolution effects on 
emission- line AGN diagnostics in §73] before summariz- 
ing our main conclusions in SjE] We assume a flat cosmol- 
ogy with f2 m = 0.3, fl\ = 0.7, and ft, = 0.7 throughout. 

2. LOW-REDSHIFT GALAXY SAMPLE 

Our low-redshift calibration s ample comes fro m the 
Sloan Digital Sky Survey (SPSS: lYoAl^aTll20Ml The 
limiting magnitude of the spectroscopic sample is r < 
17.7. Our analysis is based on data products from Data 
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Release 4 (|Adelman-McCarthv et al.l [20061) . namely the 
value-added galaxy catalogs 1 from the Max-Planck Insti- 
tute for Astronomy (Garching) and John Hopkins Uni- 
versity. In these catal ogs, the stellar masses are cal- 
culat ed as described in Kauffm ann et al.l ()2003l ) assum- 
ing a Kroupal (|2001l ) initial mass function (IMF), while 
emission line measure ments follow the procedures from 
ITremonti et~aT1 ([20041 ). 

In order to avoid strong aperture bias due to SDSS 
fiber size, we constrain the redshift range to 0.05 < z < 
0.2. This lower limit c orresponds to a minimu m cov- 
ering fraction - 30%. IKewlev fc Ellison] (|2008l) found 
that using a covering fraction down to 20% (z ~ 0.04) 
is not sufficient to obtain global emission line properties 
for massive galaxies (with M* > 10 10 M Q ). That is why 
we adopt a stricter requirement here. 

Imposing a signal-to- noise ratio cut S/N> 3 on all 
emission lines to be used in the diagnostic diagrams 
([On]A3727, Hft [Om]A5007, Ha, [Nn]A6584 and 
[Sn] AA6717,6731; g3H]), we obtain a sample of 110,205 
emission-line galaxies. 

3. EMISSION-LINE DIAGNOSTIC DIAGRAMS 

Here we aim to create a modified version of the BPT 2 
diagram, involving [Nil] A6584/Ha and [O in] A5007/H/3, 
by replacing the ratio of the redder emission lines 
QN n] A6584/Ha) because they shift into the near- 
infrared at z > 0.4. To that purpose, we introduce 
the Mass-Excitation (MEx) diagram below, which we 
will subsequently apply to intermediate redshift galax- 
ies (0.3 < z < 1) in 

3.1. Calibration Using z ~ 0.1 SDSS Galaxies 

In this section, we introduce and calibrate the Mass- 
Excitation diagram and compare it to one of the orig- 
inal BPT diagrams. The emission line ratios used in 
the best-known version of this diagram ([Nil] A6584/Ha 
and [O in] A5007/H/3) probe a combination of the ion- 
ization parameter and the gas-phase metal abundance 
within galaxies. As shown in Figure [Ha), SDSS galaxies 
form a well-defined excitation sequence on the lower-left 
of the BP T diagram (below the se mi-empirical dividing 
curve from Kauffmann et al.l [20031 solid line), while the 
galaxies containing AGN form a plume extending to the 
top right part of the diagram. The higher ionization pa- 
rameter and/or harder ionizing radiation that occur only 
in the presence of AGN cause the line ratios to lie above 
and to the right of t he maximum starbur st curve (dashed 
line) developed by K ewlev et al.l (|2001[) . Galaxies that 
are located between both curves are believed to host a 
mixture of star-formation and AGN and are sometimes 
called composites. We adopt this nomenclature in the 
remainder of this work. 

While the BPT diagram shown in Figure QJa) is used 
extensively to identify the source of ionization in galax- 
ies, the [N n] and Ha emission lines become unavailable 
in optical spectra at z > 0.4. What would be a good 
substitute for [Nil] /Ha? The [Nil] /Ha line ratio pro- 
vides an indicatio n of the gas-phase meta llicity in star- 
forming galaxies (jKewlev &: Ellison! 120081) . The empiri- 

1 http:/ /www. mpa-garching.mpg.de/SDSS/DR4/ 

2 Named after the last names of the three authors wh o intro- 
duced it: Baldwin-Phillips- Terlevich HBaldwin et al.1119811) 



cal mass-metallicity relation fe.g.. ITremonti et al.ll2004T : 
ISavaglio et al.l [20051) suggests a physical connection be- 
tween that line ratio and the stellar mass, as displayed in 
Figure [Tfb) . While the M* — Z relation applies to star- 
forming galaxies without AGNs, there is another effect 
that makes stellar mass a good choice. The [Nil] /Ha 
ratio saturates at high values for normal star-forming 
galaxies and only more extreme conditions such as those 
encoun tered in the presenc e of an AGN can yield larger 
values ([Kewlev et al.l[2006l IStasinska et al.l 120061 ). Be- 
cause AGNs tend to be found in hosts with high stellar 
mass (Kauffmann et al. 2003), these systems have both 
larger [Nil] /Ha and M* values. This feature puts them 
in a location of the parameter space of the modified di- 
agram that is analogous to their original location on the 
BPT diagram with respect to purely star-forming galax- 
ies (i.e., higher and to the right). Consequently, the new 
AGN diagnostic is obtained by substituting stellar mass 
for the redder emission line ratio. 

Indeed, we find that the BPT-SF and BPT-AGN 
classes are well separated on the new Mass-Excitation 
(MEx) diagram (Figure [lie)). We define two empiri- 
cal dividing lines that maximize the separation between 
galaxy classes, especially between the BPT-AGN class 
(above and to the right of the lines) and the purely star- 
forming BPT-SF class (below and to the left of the di- 
viding lines). We note that the location of the BPT- 
composite galaxies on the MEx diagram overlap with 
galaxies belonging to the other classes as shown in Fig- 
ure[IJd). Their locus peaks in the region between the two 
empirical curves, which we dub the MEx-intermediate re- 
gion. We note that 48% of the galaxies in this region are 
BPT-compositcs. 

The number of galaxies of each BPT-class are reported 
in Table Q] for the three regions of the MEx diagram, 
and the main empirical dividion (top curve) is defined as 
follows: 

_ ( 0.37/(x - 10.5) + 1. ifx<9.9 , . 

1 ao + a\x + a2X 2 + a^x 3 otherwise, ^ ' 

where y = log([Om] A5007/H/3) and x = log(M*). 
The coefficients are the following: {00,01,02,03} = 
{594.753,-167.074,15.6748,-0.491215}. Similarly, the 
lower curve defined as: 

y = 800.492 - 217.328x + 19.6431x 2 - 0.591349x 3 , (2) 

and is used over the range 9.9 < x < 11.2. 

Using the BPT classes as our reference, we compute the 
completeness and contamination fractions of the MEx se- 
lection for both the AGN and SF classes, separated using 
Eq. [T] The fractions are displayed in Figure [2] For the 
MEx-AGN class, we find that the comple teness is very 
high in terms of BPT- AGNs (i.e., above the lKewlev et al.l 
(|2001l) line), reaching close to 100% at high stellar mass 
(see the red diamonds in the top panel of Figure [5]). The 
fraction of BPT-composite galaxies that are selected in 
the MEx-AGN side rises more slowly with stellar mass: 
from -20% at M* < 10 10 M Q to 95% at M* > 10 115 M Q 
(black triangles). In the top panel, the contamination 
fraction is defined as the fraction of BPT-SF galaxies 
(below the iKauffmann et al.l (|2003l ) line) in the MEx- 
AGN side. The contamination fraction peaks at 10-20% 
for galaxies with M+ ~ 10 10 M , but it mostly stays well 
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Fig. 1. — Emission line di agnostic diagrams : (a) B PT diagram: purely star-forming galaxies (SF, shown in blue) form a sequence below 
the solid line adapted from Kauffmann et al. (2003) whereas galaxies hosting an AGN (red) tend to occupy the area above and to the 
right of the dashed curve IKewlev et a.1.11200111 . Galaxies located between the two curves are composites, i.e., having mixed SF and AGN 
contributions (BPT-comp., green), (b) Stellar mass as a function of the line ratio [N II] A6584/Ha. For star-forming galaxies the trend 
is analogous to the M* — Z relation (solid contours). For the total sample including composites and AGNs, there is a stronger increase 
in [Nil] /Ha compared to stellar mass but AGNs reside in galaxies with both a high value of [Nil] /Ha and a high stellar mass. These 
two features indicate that stellar mass is a viable substitute when [N II] A6584 or Ha are not available. The modified diagnostic diagram 
is shown in (c) for the galaxies classified as SF or AGN on the BPT diagram (blue and red contours, respectively) and in (d) for those 
classified as composites (green contours). The MEx-intermediate region, located between the two empirical curves on the MEx diagram, 
contains a mix of composites (BPT-comp.) and star-forming (BPT-SF) galaxies. In all panels, the contours indicate the density of points 
(in bins of 0.075 dex X 0.075 dex) and are logarithmic (0.5 dex apart, with the outermost contour set to 10 galaxies per bin). (A color 
version of this figure is available in the online journal.) 



below 10%, with an overall fraction around 6% (blue as- 
terisks). Global completeness and contamination values, 
calculated for the entire SDSS sample (all stellar masses), 
are shown with the larger plotting symbols on the right 
hand side of the figure. 

Similarly, the completeness for the MEx-SF selection 
is defined as the fraction of the BPT-SF galaxies that are 
correctly identified. The completeness is close to 100% 
(blue asterisks in Figure 0b]) over a wide range of stellar 
masses with a drop off at M* > 10 11 Mq. We note 
that at such high masses, there are very few purely star- 
forming galaxies (~200 and <10 in the last two bins) 
whereas there are of the order of 10 4 galaxies per bin at 
lower mass (Figure [5Jc] ) . The contamination with the 
BPT- AGN class is extremely low (red diamonds) , on the 
order of 0.3%. However, the contamination fraction for 
BPT-composite goes up to -30% at M* > 10 11 M Q . 
These galaxies lie mostly in the intermediate region of 
the MEx diagram. Whether they should be included 
with the star-forming sample or rejected will depend on 



the specific science goal. 

We report the numbers for the MEx-intermediate class 
(i.e., between the lines defined by Eq.Q]and[2]) separately 
in Table [TJ Even though the AGN contamination of the 
MEx-SF class appears to be more substantial at high 
stellar mass, the effect is exaggerated by the small num- 
ber of purely star-forming galaxies with M* > 10 11 Mq. 
In fact, the AGN completeness peaks at the same mass 
range where the MEx-SF class appears the most contam- 
inated, indicating that overall, the majority of the most 
massive galaxies harbor an AGN. Conversely, very few 
AGNs reside in low-M* hosts. 

The scarcity of BPT- AGNs in host g alaxies with M+ < 
10 10 Mp) was previously noted by iKauffman n et al.1 
(2003). These authors found that adding the emission 
from low-luminosity AGNs (with 10 5 < £[onr]A5007 < 
10 6 Lq) to low-M* star- forming galaxies would signifi- 
cantly alter their line ratios and move the correspond- 
ing points into the composite or AGN regions of the 
BPT diagram (in 93% of the cases with low-luminosity 



AGN Diagnostic at z > 0.3 



5 



1.0 ■ 



E 0,8 

o 

o 

<*j 0.6 



0.4 ■ 



S 0.2 - 



0.0 



S i.o- 

g" 0.8 - 

E 

o 

c« 0.6 

a 

o 

£ 0.4 h 



5 0.2 ■ 



(a) KIEx AGN 



(b) MEx SF 




O AGN ■ 



^SF/AGN. 



X SF 
'I 



* SF 



AGNs, and >99% for high-luminosity AGNs, i.e., with 



o.o 



10000 



1000 - 



100 ■ 



10 ■ 



10 11 

log M. [M a ] 

Fig. 2. — Completeness and contamination rates in bins of stellar 
mass (0.25 dex bins). Top: For the MEx-AGN selection, the AGN 
completeness is denned as the percentage of BPT-AGNs (above 
the Kewley line in the BPT diagram) that are above the MEx dia- 
gram dividing line (red diamonds). The definition is similar for the 
BPT-composite galaxies (black triangles). The contamination rate 
corresponds to the percentage of galaxies in the AGN side of the 
MEx diagram that would be purely star-forming galaxies accord- 
ing to the BPT diagram (blue asterisks). The overall values for the 
entire stellar mass range are shown on the right hand side of the 
figure. Middle For the MEx-SF selection, the completeness is de- 
fined as the percentage of BPT-SF galaxies (below the Kauffmann 
line) that are correctly identified in the star-forming side. The con- 
tamination rates are computed separately for the BPT-AGNs (red 
diamonds) and BPT composites (black triangles) and correspond 
to percentage of the number of galaxies in the star-forming side of 
the MEx diagram. Error bars are Poissonian. Bottom Distribu- 
tion of stellar masses for galaxies that are classified as star-forming 
(blue), AGN (red), or composite (black). (A color version of this 
figure is available in the online journal.) 
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So far, we have used the BPT spectral types (SF, AGN, 
composites) as references to quantify the completeness 
and contamination rates of the MEx selection. Now, we 
turn the situation around and we show a graphical com- 
parison of the MEx diagram selection (eq. 1) mapped 
back onto the usual BPT diagram. FigureEJa) shows the 
classification on the MEx diagram with MEx-AGN and 
MEx-SF galaxies colored in orange and purple, respec- 
tively. The contours show the density of points on a loga- 
rithmic scale. As can be seen in panel (b), our new AGN 
selection (orange contours) picks all the AGNs from the 
BPT diagram (above the dashed line in panel (b)) as well 
as a fair number of BPT-composites (mostly at higher 
[Om] A5007/H/3 ratio), whereas our new MEx-SF selec- 
tion (purple contours) captures the BPT star-forming se- 
quence extremely well, with an extension into the BPT- 
composite region between the two dividing lines , espe- 
cially in the region close to the [Kauffma nn et al.l (2003) 
solid line. Galaxies in the MEx-intermediate region (light 
purple dots) are distributed over the BPT-comp. and 
BPT-SF regions. 

3.2. A closer look at LINERs 

The nature of LINERs remains controversial. Histori- 
cally, they have been identified by the unusual strength 
of their na rrow optical emi ssion lines with low-ionization 
potential (Hcckman| 119801 ) . While most st udie s con- 
sider that LINER s are accretion-powered (|Hd 119991 : 
iKewlev et al.ll20C)6|) . there are also claims that, in some 
cases, the powering source could be entirel y stellar (e.g., 
from post AGB stars and/or white dwarfs, iBinette et all 
[1991 IStasihska et all [2008) or a combination of pro- 
cesses including shocks pleckmanlll980l ) . lEracleous et al.l 
(2010) argue that AGN activity does not provide enough 
energy to produce the LINER emission in half of their 
X-ray-selected sample of 35 LINE Rs. Such a deficit 
had been suggested previously (e.g.. lHo et al.lll99"3l) . and 
could be compensated by either obscuration of UV pho- 
tons from the AGN, or by contributions from alternative 
sources such as post- AGB stars. 

Chandra X-ray observations have been used to look 
for X-ray cores within LINERs. The detection frac- 
tion is typ ically high , ranging from 50% to 70% (see the 
review bv 151120081 an d refer ences t herein ) . Employ- 
ing Spitzer spectroscopy. iDudik et al.l <|2009f ) searched for 
high- ionization lines associated with AGN activity ([Nc 
V] 14 and 24/Ltm). They find that 39% of their sample of 
67 LINERs have such detections and that many of these 
cases lack AGN signatures at optical and X-ray wave- 
lengths. They also show that the optical identification 
(from, e.g., broad Ha lines) is more subject to fail at 
higher infrared luminosities, suggesting that some AGNs 
may be missed due to dust obscuration. Their AGN 
fraction in LINERs goes up to 74% after combining di- 
agnostics in all three wavelength ranges (X-ray, optical, 
mid-IR). 

The BPT diagram that includes [O in] A5007/H/3 
against [S n] AA6717, 6731/Ha, hereafter the [Sn] dia- 
gram, can be used to tell apar t the LINERs fro m the 
Seyfert 2 (Sy2) galaxies (e.g.. IKewlev et al.l 12006). We 
apply this diagram to our SDSS sample in Figure |4j 
Panel (a) shows the [S n] diagram with dividing lines be- 
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Fig. 3. — (left) Distribution of SDSS galaxies on the MEx diagram with our empirical divisions. Galaxies above the lines as classified as 
MEx-AGN (orange) whereas the galaxies below the lines are classified as star-forming (MEx-SF, in violet) and galaxies between the lines 
are MEx-intermediatc (MEx-interm., in light purple). Using the same color-coding, we show these galaxies on the standard BPT diagram 
(right). Here, galaxies that are in the MEx-intermediate region are shown with dots while the MEx-SF and MEx-AGN subsamples are 
overplotted with contours (violet and orange, respectively). The contours indicate the density of points (in bins of 0.075 dex X 0.075 dex) 
and are logarithmic (0.5 dex apart, with the outermost contour set to 10 galaxies per bin). (A color version of this figure is available in 
the online journal.) 



tween star-forming, Seyfert, and LINER populations as 
labeled. Seyfert 2 (in red) and LINER (in orange) galax- 
ies as denned in the [S ii] diagram overlap slightly on the 
[Nil] (panel b) and MEx (panels c,d) diagrams. BPT- 
composites (between the lines on the BPT [N ii] diagram) 
fall closer to the LINER than to the Sy2 distribution. 

Composite galaxies stand out more in their [Nil] /Ha 
than in their [Sii]/Hq ratio. Consequently, the star- 
forming selection (blue contours) extends into the com- 
posite and AGN regions of the BPT and MEx diagrams 
Figure lUJb) and (c)]. As mentioned previously, the 
S ii]/Ha ratio has the advantage of splittin g the LINERs 
from t he Seyfert AGNs [dashed line from IKewlev et alJ 
(2006)]. This allows us to see that the MEx diagram 
selects both of these types of AGNs, and that the SDSS 
sample used here contains more LINERs than Seyfert 2's. 

In the remainder of our analysis, we will use a scheme 
that combines the most trusted features of each diagram 
(i.e., the [Nil] and the [Sii] diagrams), which we intro- 
duce in the next section. 

3.3. A Probabilistic Approach to Galaxy Classification 

As was shown previously, some galaxy spectral classes 
overlap on the MEx diagram [Fig. [T[b,c); [3]and|l[c,d)]. 
Namely the MEx-intermediate region contains BPT-SF 
and BPT-composites, and there is also overlap between 
the BPT-LINERs and other AGN subclasses. To better 
assess the classification in such cases we present a scheme 
based on the probability of each spectral class given a 
galaxy's position on the MEx diagram. This approach 
is useful to discriminate between star-forming galaxies, 
composite galaxies, LINERs and Seyfert 2's. 

In order to quantify the probabilities, we use a low- 
rcdshift 3 SDSS sample as a calibration set. All of the 

3 0.05< z <0.1; The upper bound of the redshift range is cho- 



BPT diagnostic emission lines are available for these 
galaxies and we can place them on the MEx diagram with 
prior knowledge of their source of ionization. The source 
of ionization is assigned according to a hybrid classifi- 
cation based on both the [N ii] and [S ii] BPT diagrams 
with the following rules: 

1. Star- Forming ( SF): galaxies below the 
iKauffmann et all (|2003l ) line on the [Nil] dia- 
gram irrespectively of their class on the [S n] 
diagram. 

2. LINER: galaxies that are classified a s AGN in the 
[N ii] diagram (above the line from IKewlev et alJ 

or in the [S 11] diagram and that are below 
and to the right of the Seyfert /LI NER dividing 
line o n the [Sn] diagram (Eq. 13 of IKewlev et alJ 
(|2006h ). 

3. Seyfert 2 (Sy2): galaxies that are classified as 
AGN in the [N 11] diagram (above the line from 
IKewlev et al.ll 2001) or in the [S 11] diagram and that 
are above and to the left of the Seyfert/LINER 
dividing line on t he [Sii] diagram (Eq. 13 of 
IKewlev et all 

4. Composite (comp): composite galaxies from the 
[N ii] diag ram (between the (|Kauffmann et ail 
120031 ) and (jKewlev et all 1200 ll ) curves) that were 
not included in the LINER and Sy2 classes de- 
scribed above. 

The resulting spectroscopically-classified sample is the 
SDSS prior sample, which we assume is representative of 

sen to remain complete in all categories including LIN ERs, whose 
detection rate decreases beyond z ~ 0.1 l|Kewlcv et aLl l2006f) . 
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Fig. 4. — (a) The [Sn] diagram. SDSS g alaxies are color-code d according to their classification as SF (blue), Seyfert 2 (red) or LINER 
(orange). The dividing lines are taken from Kcwlcy ct al. (2006). All three classes are shown on the BPT diagram (b) with the same color 
scheme. The [Sn] diagram is not very sensitive to composite galaxies (BPT-comp.) as most of them are classified as SF (blue contours). 
Similarly, the three [S II] diagram classes are shown on the MEx diagram (panel (c) for the SF and Seyfert 2 galaxies, and panel (d) for the 
LINERs relative to the Seyfert 2's). The upper empirical line on the MEx diagram selects all Sy2's and LINERs on the AGN side (above 
the lines), with some overlap between the two classes. In all panels, contours indicate the density of points (in bins of 0.075 dex X 0.075 
dex) and are logarithmic (0.5 dex apart, with the outermost contour set to 10 galaxies per bin). (A color version of this figure is available 
in the online journal.) 



galaxies out to z = 1 (see £ 17.51 for a discussion of possible 
evolutionary effects). 

Relative to the previous diagrams, the combined classi- 
fication method yields a slightly greater number of AGNs 
than each of the [N n] and [S n] diagram taken separately 
(because we use the union of the two AGN subsamples). 
Correspondingly, the composite class is slightly less nu- 
merous because some composite galaxies from the [Nil] 
diagram are selected as AGN based on the [S n] diagram 
(mostly in the LINER region). Another distinction is 
that the LINERs from the [S n] diagram that extend into 
the star- forming side of the [N n] diagram are now in the 
star-forming category. The latter category is identical to 
the original [N n] star-forming selection. 

We now have a spectral classification scheme that ex- 
ploits the strengths of both the [N n] and [S n] diagrams, 
and contains four mutually exclusive categories: star- 
forming (SF), composite (comp), LINER and Seyfert 2. 
We employ the distribution of the SDSS galaxies on the 
MEx diagram with a priori knowledge of their BPT clas- 
sification to compute the fraction of galaxies of each cat- 



egory across the MEx plane. 

For galaxies at higher redshift for which only 
[Oni] A5007/H/3 and stellar mass (M±) are available, we 
use the rectangular region on the MEx diagram defined 
by the one-sigma uncertainties on these two measure- 
ments. We compute the number of SDSS prior galax- 
ies in each category (star-forming, composite, LINER, 
and Seyfert 2) and we normalize by the total number 
of SDSS prior galaxies within the rectangular box. The 
fractions are converted to percentages. Given that the 
four classes described above are mutually exclusive, the 
sum P(SF) + P(comp) + P(LINER) + P(Sy2) = 100%. 
For example, if a region surrounding given values of 
[O in] A5007/H/3 and M* (defined by the one-sigma er- 
ror bars) contains 20,000 SF galaxies; 10,000 composites; 
10,000 LINERs and no Seyfert 2s in the SDSS prior sam- 
ple, the assigned probabilities would be P(SF) = 50%, 
P(comp) = 25%, P(LINER) = 25%, and P(Sy2) = 0%. 
Thus, this probabilistic AGN classification scheme has a 
built-in uncertainty. This is a useful feature compared to 
alternative diagrams where there is often no knowledge 
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of the reliability of a certain classification. With this new 
approach, we know whether a given galaxy is near a di- 
viding line or whether it is located far into the AGN or 
star-forming locus. In this paper, we assume that com- 
posites, LINERs, and Seyfert 2's all host an active nu- 
cleus and often use: P(AGN) = 1 — P(SF) (equivalent 
to P(AGN) = P(comp) + P(LINER) + P(Sy2)). 

The empirical division on the MEx diagram intro- 
duced earlier (§ 13. ip traces well the observed transition 
between SDSS galaxies that host AGN activity of any 
category (composite, LINER, or Seyfert 2) and galaxies 
that are most likely star- forming [Figure |5Ja,b)]. Fig- 
ure [5] shows the probabilities P(SF) and P(AGN) as a 
function of position on the MEx diagram. The lower di- 
viding curve delineates the separation between a cleaner 
star-forming galaxy sample (below) and the MEx region 
where P(AGN) > 30% (above). The region between the 
two curves contains a mixed BPT-SF/BPT-compositc 
population. We adopt the term MEx-intermediate to de- 
scribe this region of the MEx diagram and the galaxies 
that are located within it. 

An alternative diagnostic diagram, developed in par- 
allel to the MEx diagram, involves the same emission- 
line ratio on the vertical axis but makes u se of rest- 
frame U — B color rather than stellar mass ( Ya n et all 
120111 hereafter Yll). We display the AGN and star- 
forming fractions on this color-excitation diagram from 
Yll in a similar fashion as for the MEx diagram [Fig- 
ure EJc,d)]. The original dividing lines are adapted from 
Yll (straight lines) and the definition is included here 
for completeness: 

log([0 in] A5007/H/3) = max{1.4 - 1.2([/ - B), -0.1}. 

In this paper, we add a curve on the diagram from Yll 
that follows the transition where the AGN probability is 
P(AGN) > 30%. The new region between this curve and 
the straight lines is analogous to the MEx-intermediate 
region of the MEx diagram and contains > 50% of BPT- 
composite galaxies. The lower curve in Figure EJc,d) is 
defined by: 

y = 12.3914 - 27.0954a: + 18.5122a; 2 - 4.02369a: 3 , (4) 

where y = log([Om] A5007/H/3) and x = (U-B) . As in 
Yll, we calculate the rest-frame U — B color, expressed 
in A B magnitudes, by using the k-correct v4_l_4 code 
from lBlanton fc Roweisl (|2007h . 

The bivariate distributions of galaxies on both the 
mass-excitation (MEx) and the color-excitation (CEx) 
diagrams are examined in more detail in Appendix IA1 

The SDSS subsample used here has been implemented 
as a reference to provide probabilities of each class as a 
function of the location on the MEx (or CEx) diagram. 
Given a stellar mass and [O in] A5007/H/3 ratio, our pub- 
licly available IDL 4 routines 5 return a probability that 
the input galaxy belongs to the SF, composite, LINER, 
and Seyfert 2 categories. The MEx diagram is well mo- 
tivated from the successes at separating and quantifying 
the overlap of BPT-classes for a large SDSS sample of 
galaxies. We will apply this new diagnostic to a sample 

4 Interactive Data Language. 

5 http://monkey.as.arizona.edu/~sjuneau/AGN-Galaxy Classi- 
fication. html 



of galaxies at 0.3 < z < 1 in the following sections. 

4. INTERMEDIATE- REDSHIFT GALAXY SAMPLE 

Our intermediate-redshift galaxy sample is based on 
observations from the Great Observatories Origins Deep 
Survey 6 (GOODS) and the All-wavelength Extended 
Groth strip International Survey 7 (AEGIS). Most of the 
analysis is performed on galaxies at 0.3 < z < 1, al- 
though we extend the range to slightly lower values when 
we use redder emission lines (such as in £15. ip . 

Optical spectra are drawn from the Tea m Keck Red- 
shift Survey 8 (TKRS iWirth et all l200l for galaxies 
in the GOODS-North (GOODS-N) field, and from the 
DEEP2 Galaxy Redsh ift Survey (hereafter DEEP2; 
iDavis et all 120031 120071 ) for galaxies in the Extended 
Groth Strip (EGS) field. These two spectroscopic sur- 
veys have similar limiting magnitudes of Rab = 24.3 
and 24.1, respectively. The former lies within GOODS-N 
(10' x 16') and includes redshifts for 1440 galaxies (1044 
galaxies with confident redshifts 0.3 < z < 1). From the 
DEEP2 survey, we only use the pointing in the Extended 
Groth Strip (centered at 14 17, +52 30) covering 120' x 
15'. There are 6,588 DEEP2 galaxies with confident red- 
shifts 0.3 < z < 1 in this pointing. 

Both sets of observa tions were obtai ned with the 
DEIMOS spectrograph (|Faber et al J 120030 at the Keck 
Observatory and reduced with the pipeline 9 developed 
by the DEEP2 team at the University of California- 
Berkeley. However, their spectral resolution and spec- 
tral range differ due to the use of different grat- 
ings (600 line mm -1 for TKRS and 1200 line mm -1 
for DEEP2). The TKRS resolution is 4A full- width- 
half-maximum (FWHM) over the wavelength range 
5500-9800A, whereas DEEP2 spectra have a resolu- 
tion of 2A FWHM with a wavelength coverage of 
6500— 9100A. These different spectral ranges mean that 
emission lines of interest are accessible over somewhat 
different redshift ranges for the two samples. H/3 and 
[Om]A5007 can be observed out to z ~1 with TKRS 
spectra but only out to z ~0.8 with DEEP2 spectra. 

A wealth of multiwavelength data is available in both 
GOODS-N and EGS. In this paper, we utilize ancillary 
Chandra X-ray data, Spitzer IRAC data, and ground- 
based optical imaging ( §4.2p . Spitzer/IRAC photometry 
is available in all four channels (available through the 
Spitzer Science Center, also see lBarmbv et ai1l2008l for 
EGS). In what follows, IRAC photometry is used to esti- 
mate the rest-frame if-band magnitude ($472]). We take 
advantage of the fact that the Ch andra X-ray coverage is 
very deep: 2 Msec in G OODS-N (lAlexander et al.ll2003h 
and 200 ksec in EGS (lLaird et alJ 120091 : iNandra et all 
2005). This provides us with an independent AGN/star- 
forming classification scheme ($5]). The sensitivity of the 
shallower data is adequate to ensure the detection of lu- 
minous X-ray sources (L2_iokcV > 10 42 erg s _1 , a nom- 
inal luminosity threshold for AGN) out to z ~ 1. Fur- 
thermore, we can detect fainter X-ray galaxies (including 
starbursts) at all redshifts in GOODS-N. 

6 http://www.stsci.edu/science/goods/ 

7 http://aegis.ucolick.org/ 

8 http:/ /tkserver. keck. hawaii.edu/tksurvey/ 

9 http:/ /astro. berkeley.edu/~cooper/deep/spec2d/ 
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Fig. 5. — Percentages of galaxies classified as star-forming or AGN (all sub-categories) as a function of location on the MEx and CEx 
diagrams (top and bottom row, respectively). The left-hand side shows the number of SF galaxies divided by total number galaxies, or 
P(SF), in each bin (0.04 dex X 0.04 dex) for the MEx (a) and CEx (c) diagrams. The AGN percentages P(AGN) are displayed on the 
right hand side, where panel (b) is the MEx diagram and panel (d) the CEx diagram developed by Yll. (U — B)o denotes the rest-frame 
U — B color expressed in AB magnitudes. In all panels, the upper lines mark the main division between AGN and star-forming galaxies, 
whereas the lower lines approximately correspond to a transition from P(AGN) < 30% (below) to P(AGN) > 30% (above). Galaxies 
between both sets of lines are a mixed SF/composite population, which we call MEx-intermediate galaxies (or CEx-intermediate in panels 
c,d). The bivariate distribution of the SDSS sample combining all galaxy classes is shown with logarithmic contours. (A color version of 
this figure is available in the online journal.) 



4.1. Emission Line Fluxes 

Emission line fluxes from the TKRS and DEEP2 spec- 
tra were measured using automated IDL routines. For 
each emission line, we fit the continuum with a linear 
relation over 30 A windows on either side of the line 
with a biweighting scheme. This technique efficiently re- 
jects outliers and is thus robust against pixels with large 
errors (e.g., due to large sky subtraction residuals). If 
the flux density at the line peak is greater than three 
times the standard deviation of the continuum, the line 
is fitted with a Gaussian function. In such cases, we 
calculate emission line fluxes in two ways. First, we inte- 
grate the resulting Gaussian fit over a wavelength range 
corresponding to ±2.5er, where a is the Gaussian width 
parameter (= FWHM/2.35). Secondly, we directly in- 
tegrate the continuum-subtracted spectra over the same 
wavelength range. 



In most cases, we use the fluxes obtained from inte- 
grating the data directly. However, some potentially 
problematic cases were flagged for visual inspection. 
Among 509 TKRS galaxies for which both the H/3 and 
[O in] A5007 emission lines passed the signal-to-noise and 
quality criteria, 61 were flagged as uncertain. Similarly, 
245 among 2,536 DEEP2 galaxies were flagged for vi- 
sual inspection. This generally occurs when the data 
are corrupted nearby the targeted emission line, or when 
the Gaussian fit is inadequate. These objects are exam- 
ined on a case-by-case basis, and the flux measurement 
is adjusted accordingly. Among a total of 306 galaxies 
that were examined, 184 were validated whereas 122 were 
flagged as uncertain and not used in subsequent analysis. 

We correct for underlying stellar absorption at H/3 and 
Ha. TKRS spectra with a median signal-to-noise ratio 
(S/N ) per pixel greater than t hree were fit individually 
using lBruzual fc Charlotl (|2003l hereafter BC03) spectral 
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synthesis models. Utilizing IDL simplefit routines (C. 
Tremonti, private communication) , we fit the continuum 
of the galaxies with a linear combination of ten represen- 
tative stellar population templates, leaving dust obscu- 
ration as a free parameter. We subtract the continuum 
to correct Balmer lines for underlying stellar absorption. 
The median corrections (± half of 84 th — 16 th percentile 
range), expressed in terms of equivalent widths (EWs), 
are 2.8(±0.9) A at H/3 and 1.4(±0.7) Aat Ha. These val- 
ues were applied to spectra that were not fit individually 
due to low signal-to-noise ratio or uncertain spectropho- 
tometry. We note that applying the median correction to 
H/3 (Ha) for spectra with an individual fit changes their 
line fluxes by 0.08 (0.03) dex r.m.s. 

In addition to using line ratios, we will use [O in] A5007 
luminosities to quantify the strength of AGN activity 
(e.g., in jJH]). For that purpose, the slit loss corrections are 
obtained by calculating synthetic photometry from the 
spectra in the band nearest to the observed wavelength of 
[O in] A5007(usually ACS F775w or F8Uw for GOODS- 
N and EGS observations, respectively) and comparing to 
the true observed photometry. The synthetic photome- 
try is obtained by applying the filter curve and integrat- 
ing over the wavelength range of interest in the observed 
band. The ratio between the total flux from observed 
photometry and the synthetic photometry is used as a 
multiplicative correction factor. The median slit loss cor- 
rection is a factor of two. 

4.2. Stellar Masses 

The stellar masses of our intermediate-redshift sam- 
ple were calculated by fitting stellar population synthe- 
sis models to spectral energy distributions (SEDs) mea- 
su red by galaxy phot ometry. The procedure is described 
in lSalim et all (|2007l) . For galaxies in EGS, the following 
photometric bandpasses are used: FUV, NU V (GALEX) 
uqriz (CFHTLS) and K (Palomar) (se e ISalim et al.1 
[20091: IGwvnll200a [20Tfl: IBundv et al.l[2006h . For galaxies 
outside of the Canada-France-Hawaii Telescope Legacy 
Survey (CFHTLS ) field-of-v i ew, w e use CFHT 12k BRI 
photometry from ICoil et all IpOOl . For GOODS-N, the 
constraints are provi ded by the fo l lowin g photometry: 
UBVRIz taken from iCapak et all (|2004l) and JK ob- 
tained with the Flamingos camera on the Mayall 4 m 
NOAO telescope. In GOODS-N, the SED fitting is 
performed for ga laxies with K < 20.5 (Vega) (as in 
IDaddi et al.H2007l but extending to lower reds hifts) . 

Th e calculations assumed a Chabrier IMF ([Chabrierl 
2003), and output a probability distribution function 
(PDF) for the stellar mass. We assume the average of 
the PDF as the stellar mass, and estimate errors from the 
2.5th and 97.5th percentiles [= (97.5PL - 2.5P£)/3.92]. 
Note that this fitting method is highly uncertain for sys- 
tems with a Type 1 AGN (identified by broad emission 
lines arising from the broad-line region) for which the 
central engine may contribute enough photons to affect 
the broad-band photometry and alter the SED fitting 
results. However, these broad-line AGNs are easily iden- 
tified and are not the object of this study, which targets 
narrow-line objects (Type 2 AGNs, LINERs, and star- 
forming galaxies). 

We note that, for SDSS galaxies, the method used 
by Salim and collaborators to derive stellar masses was 



tested against the results from lKauffma nn~et al.l (|2003f) . 
ISalim et all ((20051) found a very good agreement between 
the two distinct methods and calculated the scatter of the 
difference to be 0.11 dex (without 3cr outliers, see their 
Fig. 1), which is smaller than the typical uncertainties for 
our galaxy sample. Thus, we do not anticipate strong 
systematic differences to be associated with the meth- 
ods used to derive stellar masses for the low-redshift and 
intermediate-redshift samples that we use. 

In general, we find that the values of stellar mass cor- 
relate well with the absolute rest-frame if-band mag- 
nitudes (Mjf). The latter are obtained by applying a 
fc-correction to the observed IRAC 3.6^m photometry. 
We calibrate the relation between M* and Mk (see Ap- 
pendix [BJ for galaxies with both of these estimates in 
order to estimate a stellar mass for galaxies lacking SED 
fitting calculations. 

Starting from 2,561 galaxies with a valid stellar mass 
from SED fitting (x 2 <7) and satisfying our emission- 
line selection, we augment our sample with 251 stellar 
masses estimated from Mk- The added galaxies had 
either missing photometry or an unacceptable SED fit 
(x 2 >7). We obtain a total sample size of 2,812 galaxies 
at 0.3 < z < 1 with stellar mass estimates and valid 
[O in] A5007 and H/3 emission line fluxes. 

4.3. X-ray Luminosity and Classification 

We convert X-ray fluxes in the hard band (2 — 8 keV; 
rest-frame 2.6 — 16.0 keV at z = 0.3 — 1) to rest-frame 

2 — 10 keV luminosities, assuming a power-law spec- 
trum w ith photon index as calculated in I Alexander et ahl 
(2003) to perform the fc-correction. The fluxes are cor- 
rected for Galactic extinction but not for absorption in- 
trinsic to each galaxy. 

Our X-ray classification is based on two criteria: (i) 
^2-iokcV > 10 42 erg s _1 ; (ii) hardness ratio 10 HR > - 
0.1 (which corresponds to photon index T < 1). X-ray 
sources are classified as AGN if they satisfy at least one 
of these criteria. Otherwise, they are classified as X-ray 
starbursts. We cannot rule out that some AGNs may fail 
both the luminosity and hardness criteria due to X-ray 
absorption or intrinsically weak emission. Those objects 
are especially interesting in the framework of this study 
given the difficulty in identifying them using solely X-ray 
observations. We keep this possibility in mind as we will 
look for AGN signatures at other wavelengths besides 
X-rays. 

Some galaxies lack a detection in the hard band (2 — 
8 keV) but are detected in the full band (0.5 - 8 keV). 
In these cases, we fc-correct the full band fluxes to ob- 
tain rest-frame 2 — 10 keV luminosities assuming an in- 
dex r =1.9. These galaxies will be marked with differ- 
ent plotting symbols when using their inferred L2-10 koV 
since those values may be limits. 

Formally, our classifica t ion sc heme differs from that 
described in IBauer et al.l |2004). Those authors used 
a different X-ray luminosity threshold (Lo.5-8 keV > 

3 x 10 42 erg s _1 ) and also included criteria based on 
inferred Hydrogen column density and the presence of 
broad (> 1000 km s _1 ) or high-ionization emission lines 
in the optical spectrum. Here, we aim for a classification 

10 Hardness ratio = (H-S) /(H+S), where H and S are the number 
of X-ray counts in the hard (2-8 keV) and soft (0.5-2 keV) bands. 
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based only on X-rays, independent from optical spec- 
troscopy. Nevertheless, we note that the resulting classi- 
fication is very s imilar for galaxies t hat overlap between 
the sample from iBauer et all (|2004f ) and that presented 
here. 

We consider X-ray detections for most of our analysis 
but we also calculated X-ray upper limits for GOODS-N 
galaxies that were selected based on their [O in] A5007 
luminosity. The X-ray upper limi t s were calculated fol- 
lowing §3.4.1 of lAlexander et al.l (|2003[ ) and assuming 
r =1.9. 

5. DIAGNOSTICS AT REDSHIFT 0.3 - 1 

Now that we have calibrated the MEx diagram at 
low-redshift with SDSS galaxies, we apply it to a sam- 
ple of intermediate redshift galaxies. Our sample, de- 
scribed in §U contains 2,812 galaxies at 0.3 < z < 1 
from the GOODS-N and EGS fields with valid H/3 and 
[O in] A5007 measurements. 

5.1. Comparison with X-ray Classification 

The validity of the MEx diagram was demonstrated by 
showing a good correspondence with the BPT classifica- 
tion in 33.11 Here we show another line of support based 
on the comparison of the MEx classification with a com- 
pletely independent scheme based on X-ray observations. 
As can be seen in Figure [^a), 85% (34/40) of the X-ray 
AGNs with valid emission line measurements (S/N >3) 
are classified as MEx- AGN (26/40) or MEx-intermediate 
(8 / 40) on the MEx diagram. Thus the MEx diagram has 
a high success rate for recovering X-ray identified AGNs. 

As for the X-ray starbursts, 50% (8/16) are classified 
as SF on our new diagram, while 19% (3/ 16) are in the 
intermediate region and the remaining 31% (5/16) re- 
side in the AGN region. However, low-luminosity AGNs 
and some heavily absorbed AGNs may appear as faint 
as X-ray starbursts so this class of objects likely includes 
these different systems in addition to genuine starburst- 
ing galaxies. Indeed, the two X-ray starbursts that are 
in the AGN region of the MEx diagram but that lie at 
sufficiently low redshift to be placed on the BPT dia- 
gram (z < 0.5) stand out in the BPT-AGN region (Fig- 
ure [7]), further confirming the presence of actively accret- 
ing black holes in these galaxies. The optical diagnostics 
are thus especially useful when the X-ray signal alone 
is ambiguous (e.g., too faint to securely identify AGNs). 
Additional support for the presence of AGN in the X- 
ray starburst class is provided for at least one galaxy 
for which we have a clear detection of the [Nev] A3425 
emission line (Figure [SJ. This transition is an unambigu- 
ous tracer of AGN activity because of it's high ionization 
potential (97.1 eV). 

Next, we use the MEx classification probability method 
described in 33~3l We calculate P(AGN) by adding the 
probabilities of any AGN category (composite, LINER 
and Seyfert 2). We will show in §0 that P(AGN) > 
30% is a useful threshold to distinguish AGNs from 
purely star-forming galaxies. When X-ray AGNs, X- 
ray starbursts, and non- X-ray detections are considered 
separately, we find that 29 among 35 X-AGNs (83%) 
galaxies have P(AGN)>30%, with an average AGN clas- 
sification probability of 77%. For the X-ray starbursts, 
which all have a low hard X-ray luminosity (£2-iokcV < 
10 42 erg s" 1 ), we find that 11 among 17 X-SBs (65%) 



have P(AGN)>30% (average P(AGN)=52%). This sug- 
gests that some of them host an X-ray absorbed or X- 
ray weak AGN. Lastly, we note that some of the galaxies 
lacking an X-ray detection lie well into the AGN region 
with P(AGN) up to 100%. These X-ray faint AGN can- 
didates are studied in more detail in 35.21 Like the X-ray 
starburst class, they may in fact include weak or heavily 
absorbed AGNs. 

Not all of the X-ray detected sources have valid emis- 
sion line measurements. X-ray galaxies for which only 
one of H/3 or [O III] A5007 is measured and the other line 
has an upper limit are shown in Figure . In these 
cases, the [O ill] /H/3 line ratio is either an upper or lower 
limit. We find that most X-ray classified AGNs appear 
consistent with their optical classification given their lim- 
its on the MEx diagram. X-ray galaxies for which no 
lines could be measured due to insufficient data quality 
or intrinsic weakness of the emission lines are considered 
in 3Q1 

We note that the locus of our intermediate-redshift 
star-forming galaxies may be offset slightly from the 
SDSS contours. A more careful analysis of the selec- 
tion effects would be required to determine whether such 
a shift is real or simply results from selection biases. We 
discuss possible evolutionary effects in 37.51 

Overall, the MEx diagram is a new tool that can be 
used on galaxies with [O in] A5007 and H/3 line fluxes as 
well as stellar mass. It is thus possible to apply a robust 
AGN diagnostic to optical spectroscopic samples out to 
z ~ 1 , provided the galaxies also have photometry allow- 
ing a stellar mass estimate. An independent X-ray clas- 
sification scheme supports the validity of our diagram. 
This means that surveys with optical spectroscopy but 
lacking X-ray coverage can still benefit from a reliable 
AGN diagnostic. 

In order to complete the comparison between the op- 
tical MEx diagram AGN classification, and the X-ray 
based classification, we perform two experiments. First, 
we use X-ray stacking to search for AGN signatures in 
galaxies that are undetected in the X-ray observations 
( 35. 2p . Second, we calculate the demographics of all X- 
ray AGNs regardless of the presence of optical emission 
lines ( 35.3[) . This exercise highlights the complementar- 
ity of the optical and X-ray selection techniques. 

5.2. X-ray Stacking 

Stacking the X-ray observations allows us to probe X- 
ray emission to fainter flux levels and, in the case of 
significant detections in more than one energy band, to 
estimate the X-ray spectral slope. The X-ray spectral 
slope (or photon index T) can be used to distinguish be- 
tween different source types. While unabsorbed AGNs 
and low-mass X-ray binaries (LMXBs) exhibit a large 
index (r > 1.3 and > 1.7, respectively), high-mass X- 
ray binaries (HMXBs) can yield flatter slopes with in- 
dices down to (r > 0.5) and absorbed AGNs have still 
flatter slopes with (—1 < T < 1). These ran ges are 
illustrated on Figure 2 of lAlexander et all (|2005l see ref- 
erences therein). 

In the case of X-ray absorbed AGNs, the flattening 
is caused by the photoelectric absorption of soft X-ray 
photons by heavy atoms in neutral gas along the line 
of sight. Because softer photons are preferentially ab- 
sorbed relative to harder ones, Compton-thick AGNs 
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Fig. 6. — Application of the MEx diagnostic to identify the presence of active nuclear activity within galaxies at 0.3 < z < 1. Contours 
show the SDSS low- 2 sample (evenly spaced on a logarithmic scale), (a) Our intermediate redshift sample is superimposed (filled squares) 
and, when available, the X-ray classification is marked with larger symbols [star symbols for X-ray starbursts; small (large) filled circles for 
X-ray AGNs with Lx < 10 43 erg s _1 (Lx > 10 43 erg s — 1 )]. We also show galaxies with only one emission line detection and one upper 
limit (b). The resulting [Oin]/H/3 ratios are upper limits when only H/3 is detected (yellow downward triangles) or lower limits when only 
[Om] A5007 is detected (orange upward triangles). The X-ray classification is marked with open triangles for X-ray starburst, and solid 
black triangles for X-ray AGNs [small (large) for X-ra y AGNs with Lx < 10 43 erg s _1 (Lx > 10 43 erg s — )]. The empirical lines on the 
MEx diagram are described in the text (see § £13.11 HOI . The error bar shown in the legend represents the typical uncertainty although we 
plot individual error bars for objects whose uncertainty are significantly larger. This diagram is applicable out to z ~ 1 and yields spectral 
classes that are very consistent with the independent X-ray classification shown here. (A color version of this figure is available in the 
online journal.) 



(with Nh > 10 24 cm~ 2 ) should exhibit a flat slope with 
I 1 < 1. However there are at least two situations where 
Compton-thick AGNs instead have a steep X-ray slope at 
energies < 10 keV. First, the shape of the X-ray spectrum 
varies greatly when the X-ray emission is dominated by 
reflection rather than direct transmission. In this case, 
the slope of the spectrum does not truly reflect the hard- 
ness of the radiation and one must use other indicators 
to estimate X-ray absorption (e.g., the equivalent width 
of the Fe Ka line). Second, the photon index at energies 
0.5 — 10 keV can be significantly altered when soft emis- 
sion from starburst activity is coincident with the harder 
emission from an absorbed AGN. We keep these possi- 
bilities in mind when we interpret the results of X-ray 
stacks. 

In this paper, we perform two different X-ray stack- 
ing analyses. The first analysis is a proof-of-concept of 
the MEx diagram. In this case, the goal is to stack the 
X-ray signal for subsamples defined from the MEx clas- 
sification. A more detailed description and the results 
are provided in this section. For the second analysis (de- 
scribed below in §6.3|) . we apply additional criteria to 
target specifically X-ray absorbed AGN candidates. 

In this section, we are interested in galaxies that are 
not detected individually in the deep 2 Msec Chandra 
observations in GOODS-N but that have the required 
emission line measurements for the MEx diagram. Em- 
ploying a stellar mass cut-off at 10 10 M Q and requiring 
H/3 and [Olll]A5007 emission line detections (S/N>3), 
there are 69 such galaxies at 0.3< z <1. 

For the X-ray stacki ng we adopt a similar approach 
to lWorslev et al.1 (J2QQ5), who stacked sources and calcu- 
lated the significance of the stacked result using 10,000 
Monte Carlo trials. In all cases, we limit the X-ray stack- 



ing to sources within 6 arcmin of the Chandra aim point 
to maximize sensitivity and we stack the sources in the 
soft and hard bands. Sources that lie within a factor of 2 
of the 90% encircled energy radius of another X-ray de- 
tected source are excluded. We only consider the stacked 
signal as a significant detection when the number of de- 
tected counts in a given band exceed the background 
count rate determined from the Monte Carlo trials by 
> 3 a. 

Following the criteria outlined above there are 47 (out 
of 69) X-ray undetected galaxies in GOODS-N for which 
we can stack the Chandra data. We select subsamples 
of the X-ray undetected galaxies to stack based on the 
probabilities P(AGN) & P(SF). Here we use P(AGN) 
as the total probability for the composite, LINER and 
Seyfert subclasses. We divide the sample at P(AGN) > 
50% and also at P(AGN) > 30% in order to include the 
MEx-intermediate region. 

In the former case, we find significant detections in 
Chandra's soft band and marginal detections in the hard 
band for these P(AGN) > 50% and P(SF) > 50% 
subsamples of 25 and 22 galaxies, respectively. The X- 
ray spectral slopes are respectively flat (r ~ 0.8) and 
slightly steeper (r ~ 1.2; Table [2]). A flat X-ray spectral 
slope of r ~ 0.8 unambiguously indic ates the presence 
of absorbed AGN activity (see Fig. 2 of lAlexander et all 
[20l)5llMaiolino et al.lll998h in the P(AGN)>50% subsam- 
ple. The stacked signal of the complementary subsam- 
ple may correspond to a mixture of star-forming galaxies 
and some absorbed AGNs. However, the difference in T 
between the P(AGN) > 50% and P(SF) > 50% sub- 
samples is not statistically significant so there is a fair 
likelihood that some absorbed AGNs are present. This 
result suggests that P(AGN)>50% may not be sufficient 
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Fig. 7.— The BPT diagram, wh ere BPT-AGNs are found 
above the dashed line (adapted from Kcwlcv ct al. 2001), star- 
forming (BPT-SF) ga laxies below the solid line (adapted from 
Kauffmann ct al. 2003), and composites (BPT-comp.) between the 
two lines. Contours show the SDSS low-z sample (evenly spaced 
on a logarithmic scale). The z ~ 0.35 galaxies for which we have 
all four diagnostic lines from TKRS (0.2 < z < 0.5) and DEEP2 
(0.3 < z < 0.4) are mostly distributed along the star-forming se- 
quence, with a few galaxies extending along the AGN plume (filled 
squares) . We identify galaxies with a Chandra detection as a func- 
tion of their X-ray classification [star symbols for X-ray starbursts; 
small (large) filled circles for X-ray AGNs with Lx < 10 43 erg s — 1 
(L x > 10 43 erg s" 1 )]. While all the X-ray identified AGNs lie 
above the Kauffmann line (solid line), there are also two galaxies 
nominally classified as X-ray starbursts that lie in the region where 
AGNs are prominent. These systems also lie in the MEx-AGN re- 
gion [Figure |6j a)] and they may be X-ray weak or X-ray absorbed 
AGNs. In either case, the optical diagnostic is a useful comple- 
ment to X-ray observations alone. (A color version of this figure is 
available in the online journal.) 

to recover all of the AGNs. 

Imposing a lower cutoff at P(AGN) > 30% yields 34 
objects to stack. Stacking these sources gives significant 
detections in the soft and hard bands corresponding to a 
flat X-ray spectral slope of T ~ 0.8, again unambiguously 
indicating the presence of absorbed AGN activity. We 
then stacked the remaining 13 systems, with P(SF) > 
70%. In this case, we only find a significant detection in 
the soft band. The steep X-ray spectral slope (T ~ 1.7) is 
consistent with the X-ray emission of these galaxies being 
dominated by star formation processes. These analyses 
are summarized in Table [5] and provide good first-order 
confirmation that the combination of the MEx diagram 
and our probabilistic approach provides a comparatively 
clean selection of star-forming galaxies and AGNs. 

Our results suggest that using P(AGN) > 30% as a 
threshold for nuclear activity leads to a cleaner separa- 
tion between AGNs and star-forming galaxies than using 
P(AGN) > 50%. We examine the objects with interme- 
diate AGN likelihood separately by stacking the 12 galax- 
ies with 30% < P(AGN) < 50%. As expected, galaxies 
in this subset are likely composite systems. We calculate 
P(comp), P(LINER), P(Sy2) separately and, on aver- 
age, the composite class is 5.5 times more likely than the 
LINER and Seyfert 2 classes taken together. We find a 
flat spectral slope of T ~ 0.6, statistically indistinguish- 
able from the P(AGN)>50% subsample. This result in- 



dicates that some galaxies with 30% < P(AGN) < 50% 
host an X-ray absorbed AGNs, therefore are an impor- 
tant population to search for Compton-thick AGNs and 
to take into account for a complete census of AGNs. 

We have demonstrated that the MEx diagnostic dia- 
gram works well at intermediate redshift (0.3 < z < 1) 
given that the majority of the X-ray AGNs lie in the 
MEx-intermediate and MEx-AGN regions rather than 
the MEx-SF region. The X-ray stacking analyses pre- 
sented in this section consolidate this result. Further- 
more, there are unidentified AGNs within the X-SB class 
(nominally X-ray starburst although some sources with 
faint luminosities are more difficult to classify unam- 
biguously) and within an X-ray undetected population. 
These can be identified using the MEx diagram. 

5.3. Demographics of X-Ray Selected AGNs 

So far, we have compared the X-ray and MEx classifi- 
cation schemes for galaxies with valid emission line mea- 
surements (S/N>3 for H/3 and/or [O III] A5007). Here, we 
also consider X-ray detected sources that have sufficient 
spectral coverage for both emission lines but nonetheless 
lack detections. We use the X-ray AGN classification de- 
scribed in E|4.3l to keep only secure X-ray AGNs (i.e., with 
^2-iofeey > 10 42 erg s -1 or hardness ratio HR >-0.1). 

Combining GOODS-N and EGS subsamples, there are 
101 X-ray detected AGNs with H/3 and [Om]A5007 
within the spectroscopic wavelength coverage. We re- 
quire S/N>3 for either H/3 or [Om]A5007 in order to 
perform the MEx classification. Based on this criterion, 
33% of the X-ray selected AGNs lack optical emission line 
signatures. Performing the MEx classification for the 
remaining 68 objects with emission line measurements 
shows that 60 are in the MEx-AGN or MEx-intermediate 
region, and the remaining 8 lie in the MEx-SF region. 
To summarize the demographics, 59% are in the MEx- 
AGN or MEx-intermediate regions (48% and 11%, re- 
spectively), 8% are in the star-forming region, and 33% 
lack a classification because their emission lines are either 
too weak or corrupted. 

As a comparison, Yll used an X-ray selected sample in 
EGS that is a superset of the EGS sample studied here. 
They combined DEEP2 spectra and MMT spectra that 
were obtained specifically to target X -ray sources tha t 
were not part of the DEEP2 slit masks (|Coil et al.ll2009ft . 
For their X-ray selected sample, they found that 51% met 
their AGN class based on the color-excitation diagram, 
22% were in the star-forming region of their diagram and 
25% lacked emission lines. Yan and collaborators refer 
to the latter as X-ray bright, optically- normal galaxies 
(XBQNGs). which they study in m ore detail (also see 
iRigbv et all 120061 : iTrump et al.ll2009f ). 

A distinction between the present work and the results 
presented in Yll is their lower signal-to- noise cutoff for 
emission line measurements (S/N>2). We prefer to use a 
more conservative selection by imposing S/N>3 based on 
the visual inspection of the spectra and the large resid- 
uals often present near the lines with S/N<3. Allowing 
for a lower S/N is likely responsible for the apparent 
smaller fraction of galaxies lacking emission line mea- 
surements in Yll although the difference between the 
two analyses is not statistically significant (~ la). We 
stress that using the intermediate region - where compos- 
ite and star-forming galaxies overlap - helps to recover 
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Fig. 8. — Two example TKRS spectra with the most obvious [Nev] A3425 detections. While the top spectrum is for a galaxy that was 
also identified as AGN in the X-rays (J123608. 13+621036. 2, ,2=0.679), the bottom one is for a galaxy that was classified as X-ray starburst 
(J123645. 40+621901. 3, z=0.455). The spectra are continuum-subtracted (grey line) and smoothed with a running median (thicker black 
line). Typical emission lines are marked with vertical dotted lines and labeled. 



a larger fraction of the X-ray detected AGNs, indicat- 
ing the importance of composite systems in achieving a 
global AGN census. 

In summary, we hnd that X-ray selected AGNs cover 
the full range of optical spectroscopic classes (AGN, SF, 
no or weak emission lines). Our results are similar to 
previous work. Also, we have shown that the MEx and 
X-ray AGN selection methods remain complementary. 
Some AGNs lack X-ray detection (or are X-ray faint) and 
we found evidence for absorption among that population. 
Given the importance of X-ray absorbed AGNs, and the 
possibility to start identifying this missing population, 
we will now focus on this topic. 

6. COMPTON-THICK AGN CANDIDATES 

While soft (0.5—2 keV) X-ray emission may be pro- 
duced by star-formation activity, a large amount of hard 
(e.g., 2—10 keV) radiation is considered to be an obvi- 
ous signature of AGN activity. However, large column 
densities (N H > 10 24 cm" 2 ) can cause severe absorption 
even for hard X-rays, making it challenging to identify 
Compton-thick AGNs even with sensitive Chandra ob- 
servations such as those used in this work. 

In the event that photons with energies between 
2—10 keV are absorbed by intervening material, it is 
possible to estimate the intrinsic luminosity of an AGN 
by looking at emission originating fro m larger scales 
than the source of abso rption (e.g., IMulchaev et al.l 
11994 IBassani et al.l [T999h . Because [O in] A5007 emis- 
sion arises in the narrow-line region, which extends sev- 
eral hundred parsecs beyond the active nucleus, it is 
not subject to the same small-scale nuclear absorption 
as hard (2-10 keV) radiation. This type of compact 
absorption is expected for geometries involving an ob- 
scur ing torus as sug gested in AGN unification mod- 
els (jAntonuccil 1 1 9931 ). Thus, the ratio of X-ray (2- 
10 keV) to intrinsic (dust corrected) [Om] luminosi- 
ties can serve as a Compton- thickness paramete r T [= 
L(2-10keV)/L([OiIl] A5007)] (|Bassani et al.l H999). This 



quantity has two main caveats. 

The first caveat is the uncertain obscuration of 
[O III] A5007 by dus t in the host galaxy (ext ernal to the 
nar row line regions , Maiolino fc Rieke 1993). For exam- 
ple, lLaMassa etaTl (|2009h found a median correction of 
1.0 mag (0.5 to 2.3 mag) at [Om]A5007 based on the 
Balmer decrement using a SDSS-selected sample of 17 
Seyfert 2 galaxies that are [O ni]-luminous. Given the 
uncertain dust obscuration corrections for [O ill] , we con- 
servatively choose to use the observed luminosities in our 
analysis while keeping in mind that a typical correction of 
1.0 mag would shift the values of the Compton-thickness 
parameter down by ~0.4 dex. 

The second caveat is the possible contribution 
of stellar photo i onizat ion to the [O in] luminosity. 
IKauffmann et al] (|2003D calculated the fractional contri- 
butions to emission line luminosities in galaxies hosting 
both star formation and AGN activity. For metal-rich 
galaxies with AGN, they found that only around 7% of 
[O in] A5007 luminosity is due to star formation (the re- 
maining 93% originates from AGN-excited gas). This is 
in contrast with other lines such as H/3 and [O n] A3727 
for which a larger contribution can be expected to come 
from H II regions rather than the narrow line region pho- 
toionised by the AGN (up to 45%-70%). 



However, the stellar contribution to L 



[Om] 



be 



much more significant in low-metallicity galaxies since 
oxygen becomes one of the main coolants and changes 
in line blanketing produce a harder ionization field 
(|Kewlev et al. 2006). To alleviate this source of contami- 
nation to the [O in] luminosity, we restrict the Compton- 
thickness analysis to high stellar mass (and presumably 
high-metallicity) galaxies, with M* > 10 10 2 Mq. 

6.1. X-Ray Absorption Versus Observed AGN Power 

In order to compare the X-ray absorption to the ob- 
served AGN power, we plot the Compton-thickness pa- 
rameter T as a function of the observed (i.e., not cor- 
rected for absorption) hard X-ray luminosity (Figure [9]) . 
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iHeckman et al.1 ((2005) found that the average values 
of log(T) are 1.59±0.48 dex and 0.57±1.06 dex for Type 
1 and Type 2 AGNs, respectively. Type Is are less ab- 
sorbed and show a tighter correlation while Type 2s in- 
clude a much broader range of X-ray absorption (going 
down to smaller values of T). We reproduce this result in 
Figure Ufa) wh e re we co mbine samples of Sey ferts from 
IHeckman etaTl (|2005[ ) and lBassani et all (|1999f ) . We also 
see a clear trend between the location of nearby Seyfert 
2's o n Figure [Ma) and the ir inferred column density Nh 
from IBassani et all (|1999[ ). As expected, the Compton- 
thick galaxies, with Nh > 10 24 cm~ 2 , tend to lie at low 
values of T [log(T)<0.25 for observed X/[0m] ratios, 
and log(T)<-0.3 after correcting Lrom] f° r dust obscura- 
tion] . 

Turning our attention to the sample of 0.3 < z < 1 
galaxies, we find a similar trend in X-ray absorption ver- 
sus hard X-ray luminosity as found for the nearby galax- 
ies. There is an apparent transition at Lx ~ 10 42 erg s _1 
(vertical dotted line). X-ray AGNs above that luminos- 
ity threshold span a restricted range in their X-to- [O in] 
luminosity ratio (1 < log(T) < 3) compared to fainter 
hard X-ray sources. Galaxies that likely host an AGN 
based on the MEx diagram lie along the lower envelope 
of the distribution, consistent with their large [OlIl]/H/3 
ratios being driven by luminous [O in] emission. Galax- 
ies that are likely star forming according to the MEx 
diagram have lower X-ray luminosities on average, and 
occupy the mid-range of the Compton-thickness param- 
eter values. 

While some of the X-ray starbursts (open star symbols) 
were also classed as SF from the MEx diagram, there are 
a few objects with a conflicting spectral class. These X- 
ray-SB but MEx-AGN galaxies (open stars around red 
circles) are clustered around log(T) = 0, the nominal 
value for Compton-thick AGNs. Given the predicted off- 
set due to using uncorrected [O in] A5007 luminosities, 
and the location of known nearby Compton-thick AGNs 
on Fig. \§l a ), we will adopt log(T) < 0.25 as a crite- 
rion for identifying Compton-thick AGN candidates in 
the next sections. In cases with low values of log(T) the 
[O in] A5007 line luminosity may be a better indicator 
of the bolometric AGN luminosity. We investigate this 
next. 

6.2. X-Ray Absorption Versus Intrinsic AGN Power 

In the previous section, we considered candidate ab- 
sorbed AGNs with a detection in the X-rays, at least in 
the full band. In this section, we will expand by including 
cases with upper limits in X-rays. Limits were derived 
for observations in GOODS-N as described in H4.3I 

We again utilize the Compton-thickness parameter in 
Figure [TU] Here, we use £[om] on the horizontal axis to 
probe the intri nsic AGN power . Nea rby [O ni]-selected 
Seyfert l's from IHeckman et all (|2005l) occupy the bright 
(elevated Lrpmi) and unabsorbed (log(T)>l) region of 
Figure [TOJa). Seyfert 2's span a broad range of X-ray 
absorption and intrinsic AGN power. 

In our intermediate- redshift sample (panel b) , we note 
that star-forming galaxies seem to occupy a region con- 
tiguous with the Compton-thick candidates, which have 
log(L 2 ~iofcev)/log(£[oiii])<0.25. To discard these sys- 
tems and keep absorbed AGN candidates, we select 
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Fig. 9. — Compton-thickness parameter T = L2— i0feev/-^[Om] 
as a function of hard X-ray luminosity (L2— lOfceV") f° r ( a ) nearby 
galaxies and (b) intermediate-redshift g alaxies. In pane l (a), t he 
[O III] -selected type 1 AGN sample from Hcckinan et al. (2005) is 
shown with blue asterisks. Type 1.5 and type 2 A GNs (green as- 
terisk s and filled circles, respectively) are taken from lBassani et al.l 
(1999). The symbols are keyed to the column density Nh in 
the case of type 2 AGNs (see legend). Dust-corrected [Om] 
values are shown with open circles for galaxies with Nh > 
10 24 cm - 2 . The intermediate-redshift sample (b) is selected to 
have log(Af* [Af0])>lO.2, and is color-coded according to the prob- 
ability to host an AGN from the MEx diagram (see color bar). The 
very likely AGNs (orange and red) tend to follow the bottom enve- 
lope of the points, and reveal some Compton-thick AGN candidates 
in galaxies t hat w ere otherwise classified as X-ray starbursts (star 
symbols; see i|4.3H . Most X-ray AGNs (points lacking star symbols) 
are characterized by log(Lx/^[o 111]) > 1- Galaxies that were not 
identified on the MEx diagram (grey points) fail the emission-line 
quality flag for H/3. 

galaxies below the ix/£[om] threshold (red line) that 
also have an significant probability of hosting an AGN 
according to the MEx diagram. We adopt P(AGN)>30% 
and stack the Chandra observations to search for hidden 
AGN signal in 3Q1 

Figure [TO] shows that the Compton-thick AGN can- 
didates, with log(T) < 0.25, span a similar range in 
[O ill] A5007 luminosity as the unabsorbed AGNs, with 
log(T) > 1. At high [Om]A5007 luminosities (e.g., 
L[om] A5007 > 10 erg s _1 ) the distribution of Compton- 
thickness values appears to be bimodal, with the AGNs 
being either largely unabsorbed (log(T) > 1) or heavily 
absorbed (log(T) < 0.25). We note that this possible bi- 
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Fig. 10. — Compton-thickness parameter T = 

^2-lOAeW-kroiit] A5007 as a function of [O III] A5007 luminos- 
ity. Both panels are shown on the same scale. Nearby galaxies 
are plotted in panel (a) and the symbols are as described in 
Fig. [9] The intermediate-redshift sample is shown in panel (b). 
The diagonal lines correspond to fixed hard X-ray (2 — 10 keV) 
luminosities of 10 40 , 10 42 , 10 44 erg s 1 as labeled. We note that 
L2-iokeV = 10 42 erg s — 1 is the nominal division that we adopt 
to separate X-ray AGNs from the X-ray starbursts although there 
are some X-AGNs below that luminosity threshold (recognized 
by their hardness ratio HR > —0.1). In addition to the X-ray 
classification, the MEx AGN probabilities are indicated with 
the color coding. P(AGN) includes all AGN classes, such that 
P(AGN)+P(SF)=100%. Galaxies that were not identified on the 
MEx diagram (grey points) fail the emission-line quality flag for 
H/3. This figure shows that star-forming/starburst galaxies occupy 
a contiguous region with the absorbed AGN candidates, stressing 
the importance of the MEx diagram in selecting likely AGNs. 
Based on the observed position of the known Compton-thick 
galaxies in panel (a), we adopt log(Lx / L[o ml )=0.25 as a criterion 
for Compton-thick candidates (horizontal line). 

modality is seen in both the nearby and the intermediate- 
redshift galaxy samples. Our sample is mostly complete 
at the bright [Om] end 11 so we do not believe that this 
trend can be attributed to selection effects. However, the 
number statistics are too small to investigate it further. 

Unsurprisingly given the fairly small volume probed, 
we mostly detect AGNs with Seyfert-like luminosities 
and below [log(L[om]) < 42.5 erg s -1 ]. However, we 
note two possible AGNs with intrinsic luminosity in the 



quasar regime [with log(Lr ni]) > 42.5 erg s^ 1 ], both 
of which appear to suffer from large X-ray absorption. 
These systems may be important testbeds for evolution- 
ary scenarios where quasars are born in a deeply en- 
shrouded environment (presumably following a major 
galaxy merger) before a blowout phase where the sur- 
rounding material is ejected r evealing the optical l y and 
X-ray bright central e ngine dSanders fc Mirabell 119961 : 
iDi Matteo etaLl [20051: iHopkins et all |2005[ ). We leave 
this analysis for future work. 

In this section we have investigated the presence of X- 
ray absorbed AGNs by using both X-ray detected and 
X-ray undetected [O in] -selected objects. We find that 
the range of intrinsic luminosities probed by Lrom] are 
comparable for Compton-thick candidates and X-ray un- 
absorbed systems. The presence of X-ray absorption is 
inferred from the X-to-[Oin] luminosity ratio (T) and 
we utilize X-ray stacking analyses in the following sec- 
tion to further justify the presence of absorption in these 
systems. 



6.3. X-ray Stacking of Highly-Absorbed AGNs 

Here we search for X-ray signatures in the X-ray un- 
detected absorbed AGNs with low X-ray/ [O in] ratios 
identified in the previous section ( 36. 2j) . There are 33 
galaxies that are part of our candidate Compton-thick 
AGN selection 12 (see Fig. Hob). Sixteen are X-ray de- 
tected and are examined on a case-by-case basis below. 
Of the 17 X-ray undetected objects, 13 lie sufficiently 
close to the Chandra aim point without lying too close to 
X-ray bright sources to allow for X-ray stacking analyses. 
We employ the same X-ray stacking method described 
in §5.21 However, while in §5.21 we were concerned with 
stacking the X-ray data for galaxies selected solely from 
the MEx classification scheme, galaxies in this section are 
selected as potential Compton-thick AGNs by requiring 
log(T) < 0.25. 

Of the 16 X-ray detected sources, six are bright at 
X-ray energies (-Fo.s-sfeeV >3xl0 -16 ergs -1 cm -2 ) and 
ten are faint (F .5-skeV <3xl0~ 16 ergs -1 cm~ 2 ). The 
galaxies with bright X-ray fluxes tend to lie at lower red- 
shift and have steep spectral slopes (r « 1.3-2.0). Two 
of the six X-ray bright sources have a clear [Ne v] A3425 
detection (Fig. [8]), an unambiguous tracer of AGN, and 
reside in the Seyfert 2 region of the BPT diagram. How- 
ever, only one of them is clearly identified as an X- 
ray AGN {L x > lO^ergs" 1 ; J123608. 13+621036.2 in 
Fig. [8]). This object is strongly absorbed (r ~0.2) and is 
a potential Compton-thick AGN. The remaining five X- 
ray bright sources fail both the luminosity and hardness 
criteria to be clearly identified as X-ray AGNs. 

As for the X-ray faint galaxies, one is detected in the 
hard band and has a flat X-ray spectral slope (r « 0.7), 
indicating that it hosts a heavily absorbed AGN. The 
other nine objects are too faint to provide significant 
constraints on their X-ray spectral slopes. However, 
from stacking the Chandra data of these nine X-ray faint 
galaxies we find a comparatively steep X-ray spectral 
slope (r pa 1.4), consistent with that found for the X- 
ray bright galaxies. 
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Fig. 11. — Compton-thickncss parameter T = L2— lOfcev/^Om] 
as a functi on of the effec t ive p hoton index (T) calcul a ted as de- 
scribed in IBassani et al.l 119991 ) and Alexander ct at] H2003T ) for 
nearby and intermediatc-rcdshift galaxies, respectively. The ex- 
pected ranges of V are shown for vario us sources of X-ray emi ssion 
at the top of the figure (taken from Alexander et at 20051, and 
references therein). While Compton-thin galaxies have T > 1 as 
expected (open circles), Compton-thick AGNs can have either a flat 
or a steep index (filled circles; see text for details). We compare the 
locus of these nearby AGNs (open and filled circles) with the 0.3 < 
z < 1 systems that we identified as Compton-thick AGN candi- 
dates based on the following criteria: log(T) < 0.25, P(AGN)>30% 
and log(M* [Af©]) > 10.2. The intermediate-redshift galaxies with 
bright X-ray fluxes (red squares) span a range of photon indices. 
One galaxy has an obviously flat slope (r ~ 0.2) and the remaining 
four have steeper slopes (r > 1) but are consistent with the range 
of values sp anned by known, ne arby Compton-thick AGNs (solid 
circles, from Bassani ct at 1999). We furthermore show the results 
from X-ray stacking of candidate Compton-thick AGNs that are ei- 
ther X-ray weak (light blue hatched region) or X-ray undetected 
(orange shaded region). The height of the shaded regions illus- 
trates the range of values for the individual galaxies used in the 
stack (upper limits in the case of the non-detections). The sample 
of X-ray non-detections clearly includes absorbed AGNs. (A color 
version of this figure is available in the online journal.) 

We stacked the Chandra data for the 13 X-ray unde- 
tected Compton-thick AGN candidates, which we treat 
separately from the individual detections. While the in- 
dividually detected galaxies show steep X-ray indices, 
for the stacked objects we find a very flat photon in- 
dex of T «0.4 (see Table 2). The flatter slope found 
here strongly suggests the presence of X-ray absorbed 
AGN activity in at least a fraction of the objects that 
were stacked. The combined low X-ray photon index 
and small values of the thickness parameter make these 
systems robust Compton-thick AGN candidates. 

Except for two objects, the results from the case-by- 
case analyses of the spectral slope are less obvious. The 
low X-ray-to- [O ill] luminosity ratio and their position 
on the MEx diagram suggests that these X-ray detected 
objects host X-ray absorbed AGN activity, which may 
appear to be in conflict with the steep X-ray spectral 
slopes found for the majority of these systems. How- 
ever, the steep X-ray spectral slopes do not preclude the 
possibility that all of these objects are heavily absorbed 
or Compton-thick AGN because the soft X-ray emission 
could be dominated by either scattered nuclear emission 



or star formation (e.g. JMatt et al.ll2000D . Indeed, known 
Compton-thick AGNs span a range of spectral slopes 
from very flat to very steep (Figure [TT|) . encompassing 
the range that we observe in our intermediate-redshift 
sample. The nearby galaxies in Figure [TT] are taken from 
IBassani et al.l (|1999D and divided between Compton-thin 
and Compton-thick at Njj = 10 24 cm~ 2 (also see Fig- 
ures [5] and [TO]) . Overall, it is therefore possible that all 
of our Compton-thick AGN candidates are genuinely ab- 
sorbed by Compton-thick material (colored symbols on 
Fig. Hip ; however, some of them cannot be unambigu- 
ously identified as such based on the X-ray slope alone. 
Higher energy observations (>10 keV) may help to con- 
firm the presence of Compton-thick AGNs in galaxies 
with a steep X-ray slope. 

6.4. Linking AGN Absorption and Optical Classification 

Next, we show the X-ray absorbed AGN candidates on 
the MEx diagram in Figure [TSTa). Interestingly, their 
location differs from that of the X-ray unabsorbed sys- 
tems [Figure [TOJb)]. The latter are distributed evenly in 
the main MEx-AGN region with a possible bias toward 
high stellar mass. In contrast, the absorbed AGN can- 
didates cluster in the MEx-intermediate region with a 
few exceptions along the low-mass and outermost con- 
tour of the AGN plume. As a reminder, the MEx- 
intermediate region is populated by BPT-composites and 
BPT-SF galaxies. This suggests that the X-ray absorbed 
AGN candidates are more likely composite systems than 
Seyferts or LINERs. 

While we consider all objects (X-ray detections and 
non-detections) together, we highlight the subset of 13 
galaxies that were not individually detected but yielded 
a very flat X-ray spectral slope of T ~ 0.4 (see green 
squares on Fig. ITOI) . The flat slope makes these systems 
robust Compton-thick AGN candidates. They occupy 
the MEx-intermediate region and the lower part of the 
AGN region where different AGN sub-classes overlap. 
The average probabilistic classification is P(AGN)= 58% 
with the AGNs three times more likely to belong to the 
composite sub-class than the Seyfert 2 and LINER sub- 
classes. This result strengthens our conclusion that com- 
posite galaxies, hosting both star formation and AGN, 
are a very important population to search for highly- 
absorbed AGNs. We further discuss these findings and 
their implications in H7.3I 

7. DISCUSSION 

Several previous studies have presented AGN diagnos- 
tic diagrams as alternatives to the BPT diagram. Be- 
cause the [Nil] A6584/Ha line ratio is only available in 
optical spectra out to rcdshift ~ 0.4, it needs to be re- 
placed in order for the diagnostic to be applicable at 
higher redshift while avoiding the need to obtain near- 
infrared spectra. In this section, we first compare the 
MEx diagram with a few other AGN diagnostics from 
the literature ($7j}. We then discuss how the MEx 
diagram contributes toward a more complete census of 
AGNs ( tj7.2p . We follow with a discussion of the impor- 
tance of composite systems in this regard ( tj7.3p and we 
compare the fraction of absorbed AGNs from our work 
to results from the literature ( iJ7.4p . Finally, we men- 
tion possible evolutionary effects that are relevant to all 
nebular line diagnostic diagrams f fc|7.5p . 
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Fig. 12.— Systems with P(AGN)>30% and log(M*[Af e ])>10.2 
shown on the MEx diagram. The points are color-coded accord- 
ing to the Compton-thickness parameter T shown in Figure 1101 
Red corresponds to the most absorbed systems, with log(T)fs;-0.7 
and the color changes gradually toward violet for the largest val- 
ues of log(T) (*s2.8). Galaxies where both H/3 and [Om]A5007 
lines are robust detections (S/N > 3; filled circles) are shown in 
panels (a) and (b) depending if they are absorbed [log(T)<0.25] 
or largely unabsorbed [log(T)>0.25]. We also show galaxies with 
detected [O III] A5007 lines and H/9 upper limits (triangles in panel 
b). In these cases, the [Om] A5007/H/3 ratio is a lower limit. The 
green squares in panel (a) mark the X-ray undetected galaxies that 
yielded a flat spectral slope when stacked. Those tend to lie in 
the MEx-intermediate region indicating that some galaxies with a 
buried AGN are also undergoing an episode of active star forma- 
tion. 

7.1. Comparison with existing diagnostics 

lLamareille et al.l (|2004l ) developed an AGN diagnostic 
diagram that, like the MEx diagram, is designed to be 
applicable to optical spectra of galaxies out to z ~ 1. In 
this case, the BPT diagram was modified by replacing 
the [Nn] A6584/Ha ratio with [On] A3727/H,3. Because 
of the greater wavelength separation between [O n] A3727 
and H/3, their flux ratio is very sensitive to dust obscu- 
ration. Thus, the authors opted for an equivalent width 
ratio to mitigate against that effect. The dividing line s 
on that diagram were recently revised (Lamareillc 2010). 

The so-called blue diagram (because it includes only 
blue lines compared to Ha) has the advantage of split- 
ting the LINERs from the Seyfert 2s. However, it suffers 
from confusion between other classes of galaxies (SF, Sy2, 
and composites, see Lamareille et al 2004; Lamareille 
2010; and Appendix C). The MEx diagram has a no- 
tably cleaner separation between SF and Sy2 galaxies, is 



more sensitive to BPT-composite galaxies, and includes 
virtually all AGNs selected in the blue diagram. 

In the case of the intermediate-redshift sample used 
here, a lot of galaxies with DEEP2 spectra have an in- 
sufficient wavelength coverage and do not contain all the 
lines required for the blue diagram. As a consequence, 
the number of galaxies for which we can use this diagram 
is much smaller than for the MEx diagram. The modest 
wavelength range required for the spectra makes the lat- 
ter more versatile than the blue diagram, which requires 
AX rest > 2900 A. 

For the galaxies that have all the required observations, 
we find that the X-ray AGN selection has a much better 
agreement with the MEx diagram (Appendix C). This 
seems to be partially due to the mixing between BPT-SF 
and BPT-Sy2 galaxies in one of the regions of the blue 
diagram, and the mixing between BPT-SF and BPT- 
composites in another region. 

The appro ach most simi l ar to the MEx diagram was 
developed bv lWeiner et al. (2007|); in fact, those authors 
used datasets that overlap with the present study. They 
combined TKRS and DEEP2 spectra in order to benefit 
from both the larger wavelength coverage of TKRS and 
the larger number of galaxies in the DEEP2 survey. They 
investigated the trends between [Nil] /Ha and [Om]/H/3 
as a function of absolute ff-band magnitude (Mh), and 
found that [Nil] /Ha was sufficiently correlated with Mh 
for the latter to be a helpful parameter when only the 
[O in] and H/3 lines are available. Weiner and collabora- 
tors did not go as far as developing a diagnostic line, but 
showed a division between the red and blue galaxies on 
those plots in the sense that the red galaxies are predomi- 
nantly occupying the region where AGNs are expected to 
lie. This suggests that optical color may be another use- 
ful discriminant between star-forming and AGN galaxies 
when combined with [Q nil /H/3 . Inde ed, this approach 
was used recently bv lYan e t al. (201T|). 

Yan and collaborators designed a similar diagram as 
presented here, except for using rest- frame U — B color in 
lieu of stellar mass. Their diagnostic produces compara- 
ble results to the MEx diagram when classifying the bulk 
of SDSS galaxies into AGN or star-forming. In detail, 
they chose to avoid targeting the composite galaxies that 
lie beyond the distribution of the more obvious AGNs 
(i.e., above the Kewley line in the BPT diagram). This 
implies that the galaxies on the star-forming side suffer 
from more frequent contamination by AGNs. Thus, we 
identified the region where the composite galaxies have 
an important overlap with star-forming galaxies (region 
between Eq. 3 and Eq. 4). The motivation is twofold. 
First, it provides a means to obtain a cleaner sample 
of star-forming galaxies. Secondly, we have shown that 
composite galaxies are an important population to look 
for highly-absorbed AGNs that are co-existing with star 
formation in their host galaxies. More detailed compar- 
isons between the MEx diagram and the CEx diagram 
developed in Yll are presented in ij3.3l and Appendix IA1 

Another alternative di agram, called the DE W dia- 
gram, was developed by iStasihska et all (|2006t ). It in- 
cludes a separation based on the 4000A-break and equiv- 
alent width of [Ne in] A3869 or [O n] A3727 (whichever 
is greater). In this case, the AGN hosts tend to have 
a larger D„(4000) with respect to normal star- forming 
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galaxies. Unfortunately, the 4000A-break is located near, 
or directly over, the detector gap for most of the DEEP2 
spectra at the redshift of interest (peaking near z = 0.7), 
so we are unable to apply this diagram for the bulk of our 
sample. The authors presented evidence that the DEW 
diagram has a better correspondence with the BPT dia- 
gram compared to the blue diagram (at least prior to its 
revision in 2010). 

One notable difference between the MEx diagram and 
the original BPT diagram (or blue diagram) is that the 
former is not scale-free. The inclusion of stellar mass as a 
parameter imposes an absolute physical scale to the prob- 
lem (whereas emission line ratios do not). This aspect 
raises questions about possible redshift evolution effects 
that could systematically affect the locus of the galax- 
ies on diagnostic diagrams. We discuss such possibilities 
below (37311 . 



7.2. A More Complete Census of AGNs 

Understanding active galactic nuclei and their role in 
galaxy evolution requires a robust AGN classification 
scheme. A complete sample should include both intrinsi- 
cally weak and intrinsically bright but absorbed AGNs. 
Weak AGNs are interesting to study and learn more 
about the low accretion phase a nd whether the AGN 
unifie d model needs to be revised l[Holl200St iTrump et al.l 
2009). Absorbed AGNs are sought to explain the unre- 
solved portion of the cosmic X-ray background, and to 
quantify their contribution to mid- and far-infrared emis- 
sion seen in infrared-luminous galaxies. They are also 
of interest to test evolutionary scenarios where galaxy 
major mergers are invoked as a mechanism to form 
quasars, through a deeply-em bedded and absorbe d pha se 
([Sanders et al.lll988t lEaSiar][T999l: iHopkins et al.ll2005h . 

The MEx diagram introduced in this paper is a tool 
to uncover AGNs that are weak in X-ray observations, 
presumably due to either intrinsic weakness or X-ray ab- 
sorption. Both of these AGN phases - weak or absorbed 
- may be found in galaxies that are simultaneously un- 
dergoing episodes of star formation, potentially masking 
AGN signatures. One important feature that allows us 
to find these systems is our probabilistic approach. 

As shown in 3Q1 using P(AGN) > 30% as a thresh- 
old provides a more complete census of AGN than using 
P(AGN) > 50%. We found that both high-probability 
AGNs (P(AGN) > 50%) and medium-probability AGNs 
(30% < P{AGN) < 50%) that are not individually de- 
tected in the Chandra data show signs of X-ray absorp- 
tion when we stacked their X-ray emission (from the re- 
sulting flat X-ray spectral slope; see Tabled]). The galax- 
ies with 30 < P(AGN) < 50% are likely composite sys- 
tems (see §5.2j) . and occupy a region of the diagram that 
overlaps with other classes, making them more challeng- 
ing to identify without the probability approach. Thus, 
identifying those AGNs is a step toward a more complete 
census of active galactic nuclei. 

Our finding that the X-ray undetected galaxies iden- 
tified using the MEx diagram are likely absorbed AGNs 
prompted us to make a more systematic search for X-ray 
absorbed systems. We probed X-ray absorption using 
the Compton-thickness parameter T as described in H6.3I 
(see references therein). We again found a significant 
detection in Chandra's hard band and a still flatter spec- 
tral slope (r ~ 0.4) for galaxies that were not detected 



individually. This result supports the idea that the X- 
ray-to-[Oin] luminosity radio does select absorbed AGNs 
as expected. We also note some galaxies with a steep 
spectral slope (r > 1) and a low value of T. These sys- 
tems are also consistent with X-ray absorption given that 
they are within the range of values of T and T spanned 
by known Compton-thick AGNs that are nearby and eas- 
ier to study in more detail (see Fig. [TT]). In these cases, 
the soft X-ray emission may arise from a scattered AGN 
spectrum or from contamination by a superimposed star- 
burst component. 

To summarize the absorbed AGN samples, we first 
found signs of X-ray absorption in MEx- AGNs that were 
not detected at X-ray wavelengths. These samples were 
selected with no a priori knowledge of X-ray absorption 
and the results are tabulated in §5.21 (Table l2|). We then 
specifically targeted X-ray absorbed AGN candidates by 
imposing log(L 2 -iofeev/£jOiii] A5007) < 0.25 (Ell). This 
yielded a sample of 33 absorbed AGN candidates: six 
have bright X-ray fluxes (and mostly a steep spectral 
slope, see red squares on Fig. [TT]); ten have faint X-ray 
fluxes and the stacked signal is characterized by a slightly 
steep slope (blue shaded region on Fig. [TT]), and of the 
remaining 17 that lack detections, 13 could be stacked 
and yielded a very flat spectral slope. So we have a sam- 
ple of 33 candidates with heavy absorption of their X-ray 
emission, including a sub-sample of 15 more robust candi- 
dates (the 13 non-detections that were stacked plus the 
two individual detections that have a flat X-ray spectral 
slope). 

7.3. Notes on Composite Galaxies 

The separation between BPT-composites and BPT- 
AGNs is likely not as sharp as illustrated on the BPT 
diagram (Figure [JJa)) and may instead be a continuum 
of fractional AGN contribution to the spectral emis- 
sion lin es used in the diagn ostic. This was suggested 
by, e.g.. iKewlev et all ([20061 ). who introduced a param- 
eter to measure the distance from t he star-forming se - 
quence on the BPT diagram ( also see I Yuan et al.l 2010) . 



The AGN branches defined in IKewlev et al.l (|2006l ) start 



from the metal-rich end of the star-forming sequence 
and follow mixing sequences toward larger values of 
[Sn] AA6717,6731/Ha or [Oi]A6300/Ha depending on 
the diagram used. The idea is that there is a contin- 
uous transition from star-forming toward AGN with an 
increasing AGN contribution. 

Like composites, which host both star formation and 
AGN, galaxies classified as BPT- AGNs (above the Kew- 
ley line) may also have a significant star formation rate 
(SFR), up to ~ 10 Afoyr- 1 ([Salim et al.ll2007h . O n aver- 
age their SFRs will be less than that i n the population of 
composite galaxies (Salim et al. 2007). One consequence 
is that the emission lines of the BPT-composites are char- 
acterized by a smaller fractional contribution from AGN 
relative to star formation, which makes the line ratios 
more ambiguous. These results also imply a continuity 
(or mixing sequence) between star-formation and AGN- 
dominatcd systems. The so-called composites may be in 
the midst of a transition between the two types. 

We might expect that some X-ray absorbed AGNs 
would reside in BPT-composite galaxies because large 
amounts of gas can provide fuel for star formation and 
act as absorbing material that attenuate X-ray signa- 
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tures. In this case the X-ray absorption could come from 
the host galaxy's ISM in addition to torus-scale absorp- 
tion. In addition, if the gas is also mixed with dust, or if 
the galaxy is viewed through a dust lane, AGN emission- 
line signatures such as [O in] A5007 can be weakened 
and furthe r diluted with emission from peripheral star 
formation (iMaiolino fc Riekd U9§5t iMalkan et all 119981: 
IGoulding fc Alexander! I2009D . causing the global signal 
to exhibit composite signatures. 

The degree to which the host galaxy ISM absorbs X- 
ray and [O in] emission may vary depending on the de- 
tailed geometry of the emitting and absorbing regions. 
The X-ray emission arises on the very small scales of the 
accretion disk, while [O in] is emitted on the larger scale 
of the narrow line regions (reaching several hundred par- 
sees to kiloparsecs from the nucleus). If the absorption 
from the host galaxy ISM were approximated as a uni- 
form screen with an extent that covers both the X-ray 
and narrow line emitting regions and with a Milky Way- 
like dust-to-gas ratio, then we would expect optical AGN 
signatures (such as [O in] A5007) to be obscured more 
than the X-ray emission. An extinction of Ay — 10 mag 
would fully obscure [O in] emission from the narrow line 
regions, but would correspond to a gas column density 
of Njj = 1.6 x 10 22 cm -2 , causing only modest absorp- 
tion of the X-ray emission. This would suppress the 
^[Om]/-^2-iOAeV ratio, leading to the opposite effect of 
the preferential X-ray absorption that we observe in some 
systems. 

However, if the foreground ISM is clumpy, then it is 
possible for a dense cloud with large Njj to obscure the 
compact X-ray emitting region, while the average ab- 
sorption to the larger [O ni]-emitting narrow line regions 
could be much less. There will be lines of sight to the ex- 
tended narrow line regions with more modest obscuration 
which will dominate the global emission line signal (rela- 
tive to more heavily obscured lines of sight). As a result, 
we may in this case expect a preferential absorption of 
the X-ray emission relative to the optical ([Om] A5007) 
emission. 

The preferential extinction of the small-scale X-ray 
emission when compared to the larger scale [O in] emis- 
sion could also be caused by a physical connection be- 
tween the gas at galaxy scales and the gas on the scales 
of the AGN torus. For ex ample, in the simulations of 
iHopkins fc Quata"er t (2010), an overall larger gas frac- 
tion on kpc scale can result in more gas funneling to the 
inner sub-pc region of the systems. These simulations 
predict that the transport of gas toward the central re- 
gion depends most strongly on the disk-to-bulge ratio 
and on the gas fractions on 100 — 300 pc scales. We note 
that these authors did not specifically model a torus. 
However, the presence of instabilities on small scales in- 
creases the gas flow to the inner regions and suggests that 
the torus, often thought to be an extension of the accre- 
tion disk, could contain more obscuring material (or have 
a larger filling factor). If this scenario were true, it would 
imply that a higher gas fraction on galaxy scales could 
lead to both a larger SFR and a more gas-rich accretion 
disk (and potentially torus), thus creating preferential 
absorption of small scale emission (hard X-rays) relative 
to larger scale emission ([O in] A5007). 

Based on these scenarios, and on our observations of 
a hard X-ray signal by stacking only 12 likely composite 



galaxies (with 30 < P(AGN) < 50%, SgT3D that were 
not detected with very sensitive Chandra observations, 
we propose that not all BPT-composite galaxies have 
intrinsically weak transition AGNs. Instead, a fraction 
of them have powerful but absorbed AGNs whose light 
is diluted with that of their host galaxies (regardless of 
whether the host galaxy ISM provides additional AGN 
absorption or not). 

In addition, the composite galaxy populations may dif- 
fer with redshift. In general, star- forming galaxies had a 
higher specific st ar formation rate (SSFR _ SFR/M+ ) 
at earlier times (jNoeske et all 120071 : lElbaz et a.1.1 12007f t . 
The larger amount of star formation at higher redhsift 
is also linked to larger reservoirs of mo lecular gas (e.g., 
iDaddi et all 120101: iTacconi et all 120101 ) . As mentioned 
before, large amounts of gas contribute to increase the 
column density along the line of sight thus X-ray absorp- 
tion. The more sizable gas reservoirs in isolated galaxies 
at higher redshift may be a way to obtain more absorbed 
AGNs without major mergers. 

There is also evidence for a larger i nfrared-to- X-ray 
lumin osity ratio with increasing redshift ([Mullanev et all 
l2010f) . further suggesting a larger SSFR in the host galax- 
ies of higher-redshift AGNs relative to the current epoch. 
The enhanced star formation may affect the line ratios 
in the sense that more galaxies will be classified as BPT- 
composites at higher redshift. While we do not study 
star formation rates in this work, we remind the reader 
that the composite galaxy population plays an important 
role in the search for X-ray absorbed AGNs and as such 
should be identifiable at high- redshift. Thus, the MEx 
diagram and the probabilistic classification scheme intro- 
duced in this Paper are expected to fulfill this need and 
to contribute to a significant improvement in identifying 
the population of missing Compton-thick AGNs. 

7.4. Absorbed AGN Fractions 

Next, we examine the absorbed AGN fraction in 
terms of the number of galaxies identified as well 
as fractional contribution to Lroml- Galaxies with 
log(Lx/i[Om])>l are considered as unabsorbed while 
galaxies with \og(Lx / L[ m]) < 0.25 are considered ab- 
sorbed and likely Compton-thick and galaxies in be- 
tween are likely absorbed but Compton-thin with 22 < 
log(A^[cm- 2 ]) < 24 [Figure [Ha)]. 

For galaxies in EGS, we compute X-ray upper limits 
using sensitivity and exposure time maps 13 . We convert 
the count rates in the ha rd band to fluxes a ssuming the 
conversion factor used bv lLaird et al.l (|2009f ) . The X-ray 
detection limit is highly variable across the EGS field. 
For the [O ni]-selected samples that we consider here, 
it varies from 6.6 x 10 -16 to 1.2 x 10 _14 ergs -1 cm -2 . 
The average value is 2.5 x 10 _15 ergs _1 , similar to the 
flux limit at which the survey is complete over 50 % of 
the area (2 x 10 -15 ergs" 1 cm" 2 : lLaird et alJl2009f) . We 
discard the X-ray upper limits of 10% (7/70) of [Oni]- 



selected galaxies with L 



O in] A5007 



> 10 ergs because 



they lie in very shallow regions of the X-ray data (with 
upper limits that are > 5 x 10~ 15 ergs _1 cm -2 ). 

However, we note that a typical upper limit value of 
2.5 x 10~ 15 ergs -1 cm -2 does not constrain the absorp- 

13 Data products from http: / /astro. ic.ac.uk/content/chandra-data- products 
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tion very tightly, especially at log(L[om] [erg s 1 ])<41.5 
where a lot of upper limits correspond to log(T)>l. 
Thus, we derive lower and upper limits to the absorbed 
AGN fraction by assuming that these galaxies are respec- 
tively all unabsorbed (i.e., they actually lie at log(T)>l) 
or all absorbed (i.e., they lie at log(T)<l). We further- 
more weigh the galaxies with the probability that they 
host an AGN according to the MEx diagnostic. The 
fractional number and [O in] luminosity contribution of 
absorbed AGNs are respectively defined as follows: 



fa 



bsorbed " 



■^N abtor 



P(AGN) Z 



Y,flT l P{AGN) 3 
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/([Oni])„ 
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where P(AGN) is the fractional probability of hosting 
an AGN varying from to 1, N a i, S orbed is the num- 
ber of X-ray absorbed AGNs, and N to tai is the to- 
tal number of AGNs (absorbed and unabsorbed). At 
log (Lfo mi [erg s _1 ])>41, the absorbed (Compton-thick) 
fractions are poorly constrained: we find /absorbed — 
25 - 81% {f Compton-thick = 12 - 81%). Using equa- 
tion[5J the fractional contribution of absorbed (Compton- 
thick) AGNs to the [O m] luminosity of all AGNs with 
log(£[om][erg s" 1 ])>41 is /([0iii]) o jw&ed =54-82% 
(31-82%). 

The constraints are slightly better with a higher [O in] 
luminosity threshold. Restricting our analysis to the 34 
galaxies at 0.5 < z < 0.8 with log(L[o ra ] [erg s -1 ])>41.5, 
we find an absorbed (Compton-thick) fraction ranging 
from 45% to 68% (17% to 68%) in number of galaxies, 
and a contribution to L[om] between 70—81% (39—81%). 
Similarly, Yll estimated an absorbed AGN fraction by 
calculating the X-ray detection probability assuming a 
constant intrinsic X-ray/ [O in] A5007 ratio and measured 
[O in] luminosities, and by comparing the expected num- 
ber of X-ray detections to the actual number of detec- 
tions at the EGS depth. Attributing the non-detections 
to X-ray absorption, they find an X-ray absorbed frac- 
tion of ~50— 60% at the same [O in] luminosity threshold 
used here, in agreement with our estimate for EGS and 
GOODS-N combined (45 - 68%). 

Using only GOODS-N galaxies provides stronger con- 
straints on X-ray absorption given the high sensitiv- 
ity of the Chandra observations in that field, but is 
subject to smaller number statistics. In this case, re- 
stricting our analysis to the 29 galaxies at 0.3 < z < 
0.8 with log(L[om][erg s~ 1 ])>40.5 yields an absorbed 
(Compton-thick) fraction of 75±18% (54-64%). For 
log(i [0l „][erg s~ 1 ])>40.5, Yll calculated that 70.5% of 
AGN fail to be detected at the EGS depth, suggesting an 
absorbed fraction (70.5 ±4.1%) that agrees with our cal- 
culations for GOOD S-N galaxies (75%). Our results ar e 
also consistent with lAkvlas fc Georgantopoulosl (2009), 
who studied 38 nearby (<70 Mpc) Seyferts and found 
an absorbed fraction 55±12%. Considering only sources 
with L 2 _ wkeV > 10 41 erg s" 1 , they find that 75±19% of 
21 galaxies are absorbed (N H > 10 22 cur 2 ) and 15-20% 
are Compton-thick. 

At quasar-like luminosities, Vigna iTet al.l (|2010[ ) found 
that 68% of SDSS Type 2 quasars at 0.3 < z < 0.8 are 



Compton-thick according to [O in] versus X-ray criteria. 
As they discuss, an [O in] -selection may bias the numbers 
toward a larger absorption fraction. Our estimates do not 
allow us to predict whether the absorbed AGN fraction 
rises or declines with [O in] luminosity. Larger samples 
with very sensitive X-ray observations would help to ad- 
dress that question. 

7.5. Possible evolutionary effects 

As noted in Sj5j there is some evidence that the 
intermediate-redshift galaxies may be offset in the MEx 
diagram relative to the lower- redshift SDSS sample. 
Other authors have reported similar shifts in the BPT 
diagram, i.e., higher- z galaxies appear displaced toward 
larger values of [Oiii]/H/3 and/or [Nn]/Ha. Some at- 
tributed this trend to varying HII region conditions 
dLiu et all 120081: iBrinchmann et all 120081: lHainline et al.l 
2009) while others claim that addi tional AGN contri- 
bution may be th e driving factor ([Groves et al.l 120061 : 
IWright et al.l[2010h. 

In particular. IWright et al.l (|2010t) presented a detailed 
study of one BPT-composite galaxy at z ~ 1.6. Us- 
ing adaptive optics combined with integrated field spec- 
troscopy, they were able to separate the AGN from the 
host galaxy emission and found that the central region 
(inner 0."2 x 0."2) occupies the AGN part of the BPT 
diagram, while the integrated measurements place this 
galaxy in the BPT-composite region and the host galaxy 
alone shares the locus of the normal star-forming galaxies 
(i.e., not shifted from the low-z sequence). This suggests 
that diluted AGN contribution may explain the offset in 
at least some cases. 

What would be the consequences of such offsets on the 
MEx diagram? If more galaxies have diluted AGN con- 
tributions, they will be selected as MEx-AGNs because 
the offset would tend to move the galaxies over the di- 
viding line. A fraction of the higher-mass galaxies are 
already subject to be BPT-composite (~ 20 — 30%) so 
in this case, moving them over the AGN line would ac- 
tually increase the completeness. On the other hand, if 
galaxies had different physical conditions in their HII re- 
gions in the past, in the sense of having larger [O m]/H/3 
ratios compared to galaxies today, then some purely star- 
forming galaxies could be moved over the AGN dividing 
line. However, the differences out to z ~ 1 are proba- 
bly small since the empirical offset is around 0.2 dex in 
log([Om]/H/3). 

Given that the situation is still under debate, and that 
selection effects are not well constrained, we do not im- 
plement an offset at this point. However, we note that 
the very good agreement between the X-ray and the MEx 
selection even out to redshift ~ 1 supports the conclusion 
that the empirical division works well over the redshift 
range considered here. Future work would be required to 
test whether it holds to even higher redshift. 

8. SUMMARY 

In this paper, we provide a new tool to gather a more 
complete census of actively accreting black holes in galax- 
ies. We argue that intrinsically weak AGNs, as well as 
heavily absorbed AGNs, are important to understand the 
connection between supermassive black holes and their 
host galaxies. We successfully find systems that belong 
to these observationally-challenging AGN classes. 
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Using 110,205 emission-line galaxies from the SDSS, 
we calibrate a new excitation diagram to identify 
galaxy nuclear activity. By combining the observed 
[O in] A5007/H/3 emission line ratio with stellar mass, 
we obtain a diagnostic applicable to galaxies out to 
z ~ 1. Here we summarize a few properties of the Mass- 
Excitation (MEx) diagram: 

1. The simplest version splits galaxies into three 
classes: MEx-AGN, MEx-SF and MEx- 
intermediate. The latter is found in a region 
where BPT-composite galaxies have significant 
overlap with star-forming galaxies. The classifica- 
tion has an excellent correspondence to that from 
the BPT diagram based on [O in] A5007/H/3 and 
[Nil] A6584/Ha. 

2. The classification scheme is refined using a novel 
probabilistic approach to predict the location of 
high-redshift galaxies on standard BPT diagrams 
based on their [O ill] A5007/H/3 ratio and stellar 
mass ( tj3.3[) . Utilizing SDSS priors, we calculate 
the probability of four mutually exclusive spec- 
tral classes (star-forming, composite, LINER, or 
Seyfert 2) as a function of position on the MEx 
diagram for each galaxy. 

3. The MEx diagram offers a more complete AGN se- 
lection than alter natives such as the blue diagram 
(lLama rcillc 2010). It is comparable to the Color- 
Excita tion (CEx) diagram developed bv lYan et al.l 
(2011) but may be particularly useful for dusty, 
strongly reddened galaxies which may have unusual 
colors. Furthermore, we add a dividing line to the 
CEx diagram to identify composite galaxies, an im- 
portant population to search for absorbed AGNs. 

4. We provide publicly available IDL routines to cal- 
culate the probabilities of different galaxy spectral 
classes based on the MEx and CEx diagrams. 

Another outcome of a nearly complete AGN selection is 
that the complementary selection results in a very clean 
star-forming galaxy sample. Indeed, with our probability 
scheme, one can weight galaxies as a function of P(SF) 
when computing average properties such as the stellar 
mass-metallicity relation or the global star formation rate 
and thus mitigate against AGN contamination. We an- 
ticipate that a wide variety of applications will benefit 
from our approach to galaxy spectral classification. 

We successfully classify AGN and star forming galaxies 
at 0.3 < z < 1 using our new diagnostic. The sample is 
drawn from the AEGIS and GOODS-N surveys. Our 
main findings are as follows: 

• An independent X-ray classification scheme shows 
that the MEx diagram selects around 82% of the 
X-ray AGNs with detected emission lines (S/N>3 
for H/3 and [Olll] A5007). When considering all X- 
ray AGNs, we find that 59% are in the MEx-AGN 
or MEx-intermediate regions, 8% are in the star- 
forming region, and 33% lack a classification (they 
fail the S/N > 3 requirement for valid emission line 
detections) . 



• The MEx diagram reveals X-ray absorbed AGN 
candidates, which were missed or mis-identified in 
X-rays while having robust signatures in their op- 
tical spectra. We combine support from three lines 
of evidence: 

1 . The subset of these galaxies with a low enough 
redshift to be placed on the BPT diagram lie 
in the BPT-composite or BPT-AGN regions 
(Fig. |7J) . In addition, one galaxy has a clear 
detection of [Ne v] A3425 in emission, a secure 
AGN tracer ({jJTTj. 

2. Stacking Chandra data of galaxies that were 
not individually detected in GOODS-N re- 
veals a hard signal and a flat X-ray spectral 
slope for galaxies that have P(AGN) > 30%. 
This suggests that at least some of the 32 
galaxies that were stacked do host an actively 
accreting black hole f H5.2p . 

3. By combining [O in] A5007 and hard X-ray lu- 
minosities, we calculate a Compton-thickness 
parameter (T = L 2 -i keV / L[ 0m ] A5007)- We 
identify 33 highly-absorbed AGN candidates. 
The presence of AGN in at least some of the 
candidates is supported by the hard X-ray 
signal (r ~ 0.4) obtained from stacking 13 
galaxies that were not detected individually 
( ij6.3p . We note that some X-ray galaxies 
have a steep spectral slope despite a low value 
of T. These objects are also consistent with 
the range of values expected from a compari- 
son with known, nearby Compton-thick AGNs 

(Fig. ED. 

• Many absorbed AGNs in the intermediate redshift 
sample are composite systems, with both star for- 
mation and active black hole accretion. The in- 
crease in the SF:AGN ratio - and hence a decrease 
in [Om] A5007/H/3 - may be related to the global 
rise in the specific star formation rate with red- 
shift. However, the details of such a connection 
remain open questions (does the host galaxy ISM 
provide additional AGN absorption? is there a link 
between host galaxy and torus properties?). The 
identification of these composites poses a consider- 
able challenge on the observational point-of-view. 
The MEx probabilistic approach proved to be key 
in this respect. 

On-going and future near-infrared multi-object spec- 
troscopy surveys will shed more light by allowing us to 
observe all the traditional emission lines directly. The 
MEx diagram will nevertheless be useful to fill in gaps 
where Ha and/or [Nil] A6584 fall in especially noisy re- 
gions of the spectra or where [N n] A6584 fails to be de- 
tected (as it is typically fainter than Ha). In principle, 
the MEx diagram could be applied out to z ~ 4 using 
near-infrared spectroscopy, thus pushing the limits be- 
yond what will be achievable when relying on Ha and 
[N n] A6584 detections (since the latter redshift out of 
the if -band at z ~ 2.5). 
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APPENDIX 

COMPARING THE MEX AND CEX DIAGRAMS 

The mass-excitation (MEx) and color-excitation (CEx) diagrams are two useful alternatives when the emission lines 
typically used for AGN diagnostics such as the BPT diagram are not available. The former is developed in this paper 
( fl3.1j) and involves using stellar mass as a discriminant between predominantly star-forming galaxies and those hosting 
an AGN. The CEx diagram was developed by Yan et al. (submitted) and involves using rest-frame U — B color rather 
than stellar mass. Galaxies hosting an AGN tend to be both massive and redder in U — B so it is not surprising that 
both these approaches give similar results for many galaxies. 

We compared the fraction of galaxies hosting AGN to the purely star-forming galaxies according to the BPT diagram 
in §3.31 We defined an additional line on the CEx diagram to provide a means to identify a region where star-forming 
galaxies are mixed with composite galaxies, analogous to our method with the MEx diagram. Here we look at the 
bivariate distributions of galaxies split into the four categories defined in §3.31 Figure [13] shows that both diagrams are 
generally comparable. The main dividing lines (upper lines) make a division between the LINERs and Seyfert 2's on 
the upper side and the star-forming galaxies on the lower side. By design, the MEx diagram selects a larger fraction 
of the composites in the AGN side. This decision was motivated by our goal to provide an increased completeness 
of AGN candidates including intrinsically weak as well as absorbed systems that may be undergoing star formation 
simultaneously with an AGN phase. We have shown the strong potential to find X-ray absorbed candidates among 
such a population of galaxies (§ §5.2|6|) . 

On the CEx diagram, we note a few galaxies in the star-forming category that seem to be outliers with very red 
U — B colors (also see Figure [3]). These may be especially dusty. Their number is small relative to the bulk of the 
star-forming population in the SDSS sample but they may be interesting targets in a different context. For instance, 
we may expect more star formation and active galactic n uclei to take place in dusty, infrared- luminous galaxies at 
higher redshift (|Le Floc'h et al.l[2005t iMagnelli et al.ll2009t ). 



24 Juneau et al. 



1.5 




Fig. 13. — Bivariate distributions of the SDSS galaxies on the MEx diagram (top row) and CEx diagram developed by Yll (bottom 
row). Galaxies are plotted separately for each classification: (a) star-forming, (b) composite, (c) LINER, and (d) Seyfert 2. The same order 
is followed in panels (e)-(h) for the CEx diagram. The lower dividing lines were added to mark the region with significant overlap between 
BPT star-forming and BPT-composite galaxies (the fraction of composites is greater than 40—50% between the lines). (A color version of 
this figure is available in the online journal.) 
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STELLAR MASS ESTIMATION 

There are 3174 galaxies with both a stellar mass estimate from SED fitting and absolute rest-frame K-b&nd magni- 
tude from observed IRAC photometry. Here, we use EGS galaxies with a SED fit to the FUV, NUV, ugriz, and K 
photometry and GOODS-N galaxies with a SED fit to UBVRIzJK bands (JO). 

The relation between stellar masses derived from SED fitting and rest-frame K-band absolute magnitudes is displayed 
in Figure [T4l The broken power-law relation is expressed as: 

logM* = -0.398 xM K + 1.357; for Mk > —20.5 AB (Bl) 
= -0.519 x M K - 1.128; for M K < -20.5 AB (B2) 

Writing the relation as M* oc L^, the slopes found in Eq. Bl & B2 imply power-law indices a = 0.99 at Mk > 
—20.5 AB and a — 1.3 at Mk < —20.5 AB. This relation is linear at the faint end but slightly steeper at the bright 
end. The steepening of the slope is presumably due to the fact that more massive galaxies are redder owing to the 
older average age of their stellar populations, and thus have higher mass-to-light ratios. 

The residuals of the broken power-law fit (bottom of Figure [T4l) have a constant scatter with luminosity. The 
dispersion of the overall sample is 0.18 dex, around the mean (median) of 0.004 (0.009) dcx. 

Errors on the stellar masses calculated using Eq. Bl & B2 are estimated by combining the average error on log(M*) 
from galaxies with SED fitting (0.12 dex) and the dispersion on the log(M+) —Mr residuals (0.18 dex). Adding these 
two contributions in quadrature yields 0.22 dex, which we use for all galaxies lacking a SED fit. 
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Fig. 14. — [Top] Stellar mass obtained from SED fitting as a function of absolute rest-frame A"-band magnitude (Mk, in AB magnitudes). 
We fit two mass-luminosity power laws for the bright and faint ends, which intersect at Mk = —20.5 AB (green lines). The dashed line 
is the extension of the fit at the faint end, with a linear slope between mass and luminosity. For reference, the yellow diamonds show the 
median log(M*) in Mk bins (every 51 points). [Bottom] Residuals of the top panel as a function of Mk- Yellow diamonds and error bars 
show the median and the 16th to 84th percentile range. There is no obvious trend in the residuals (median=0.009 dex), and the overall 
dispersion (f» (84PL — 16PL)/2) is 0.18 dex. (A color version of this figure is available in the online journal.) 
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COMPARISON WITH THE BLUE DIAGRAM 

The blue diagram discussed in §7.11 was developed in lLamareille et al.l (j2004T ) and recently improved in iLamareillel 
(2010). The diagnostic employs lines at blue rest-frame wavelengths, from [On]A3727 to [Oiii]A5007, in order to 
facilitate its use out to z ~ 1. This motivation is similar to that which guides the design of our MEx diagram, 
although our results differ significantly. 

In this diagnostic, the abscissa is an equivalent width ratio because the spectral separation between the [O n] A3727 
and H/3 lines is significant and their line flux ratio would thus be very sensitive to dust attenuation. Note the inability 
to distinguish between SF and Seyfert 2 galaxies in the region marked as SF/Sy2. We show that this ambiguity is 
completely removed with our new diagram. 
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Fig. 15. — Distribution of SDSS galaxies on the blue diagram: [O Hi] A 5007/H/3 line flux ratio against [On] A3727/H/3 equivalent width 
ratio. This diagram and dividing lines are adapted from Lamareillc (2010). Most star-forming galaxies (SF) occupy the bottom left of the 
plot (blue dots) but include a number of composite galaxies. The black dots correspond to AGNs (Seyfert 2s and LINERs), and the light 
blue points mark the galaxies in the region where SF and Sy2 classes overlap. These classifications are compared to other diagnostics: (b) 
[Nil] BPT diagram, (c) [Sir] BPT diagram, and (d) MEx diagram. The blue diagram star-forming galaxies include a number of composite 
galaxies as shown on the BPT diagram (b), as well as AGNs that are mostly LINERs according to panel (c). The ambiguous SF/Sy2 
galaxies (light blue) are well separated on the other panels, especially (b) and (d). In all panels, the contours indicate the density of points 
(in bins of 0.075 dex X 0.075 dex) and are logarithmic (0.5 dex apart, with the outermost contour set to 10 galaxies per bin). (A color 
version of this figure is available in the online journal.) 



We examine the selection functions built in the blue diagram against the three other diagnostics introduced earlier. 
On Figure [T5l the blue points represent the star-for ming galaxies (pa nel (a); note that there is a known overlap 
with the composite galaxy population as discussed in Lamareillc 2010). The pale blue points mark the ambiguous 
region where SF and Seyfert 2 galaxies are indistinguishable (SF/Sy2) whereas the black dots indicate the AGNs 
(encompassing both Seyfert 2's and LINERs). Panels (b)-(d) show the location of the same galaxies over the other 
diagnostic diagrams. The main features are that (i) the ambiguous SF/Sy2 region (pale blue points) is resolved in the 
other diagrams, which all break the degeneracy observed in the [On] A3727/H/3 EW ratio of low-M* SF galaxies and 
Seyfert 2s; (ii) the blue diagram AGNs (black dots and black contours) include some BPT-composites: more so than 
the [S n]-diagram shown in (c), but less so than the BPT (b) and MEx (d) diagrams; and (iii) a large fraction of SF 
galaxies on the blue diagram are in fact composites or AGNs according to the BPT diagram (b). 

The MEx diagram (Figure [6]) is applicable to 2,812 galaxies out to z ~ 1, whereas the emission lines in the blue 
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diagram (FigureHU) are only available for 826 galaxies (29%). While 423 out of 531 (80%) GOODS-N galaxies with H/3 
and [O in] A5007 also have a valid measurement of the [O n] A3727 line, the situation is very different for EGS galaxies. 
The DEEP2 spectra used for these galaxies span a more restricted range in wavelength and so only 403 among 2,536 
objects (16%) with valid [Om] A5007 and H/3 also have a valid [On] A3727 line flux. 

The star-forming region of the blue diagram (Figure ITS"]) accounts for 10/13 (77%) of the X-ray starbursts but also 
for 6/11 (55%) of the X-ray AGNs. The Seyfert 2 region only accounts for 3/11 (27%) of the X-ray AGNs. The smaller 
number of galaxies on the blue diagram hinder a quantitative comparison with the MEx diagram. Nevertheless, it 
appears that the former is less reliable at selecting X-ray AGNs. 
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Fig. 16. — The blue diagram, where AGNs are found above the solid curve and further separated between Seyfert 2 and LINER classes, as 
labeled. Star-forming (SF) galaxies lie below the curve including a region where they are indistinguishable from Serfert 2 galaxies (labeled 
as SF/Sy2). In principle, this diagram is applicable out to % ~ 1. However, the DEEP2 spectra used in this work span a narrow wavelength 
range and very few cover all three lines required in this diagram. While performing well at selecting the X-ray starbursts (star symbols), 
this diagram misses the majority of X-ray AGNs (filled black circles) including the brighter ones (larger circles). Contours show the SDSS 
low-z sample (evenly spaced on a logarithmic scale). Our intermediate redshift sample is superimposed (orange points) and, when available, 
the X-ray classification is marked with larger symbols [star symbols for X-ray starbursts; small (large) filled circles for X-ray AGNs with 
Lx < 10 43 erg s — 1 (Lx > 10 43 erg s )]. (A color version of this figure is available in the online journal.) 
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TABLE 1 

Demographics of the MEx Diagram 



BPT type 


MEx-AGN 


(96) 


MEx-Interm. a 


(96) 


MEx-SF 


(%) 


BPT-SF 


1465 


6.0 


9153 


51.7 


64243 


94.4 


BPT-compositc 


9782 


40.0 


8468 


47.8 


3760 


5.5 


BPT-AGN 


13212 


54.0 


90 


0.5 


31 


0.04 


All 


24459 


100 


17711 


100 


68034 


100 



a MEx-intermediate region: between the two curves defined by Equations [T1 and 

13 



TABLE 2 

X-ray Stacking of Non-Detections 



Sample 


N a 


< z > 


Soft b 


Soft 
a 


Hard b 


Hard 
a 


HR C 


r d 


Comment 


(i) MEx Diagram Selection (see 95.21 


) 














P(AGN)>50% 


22 


0.70 ± 0.18 


1.61 ± 0.54 


4.6 


1.61 ± 0.81 


2.9 


0.00 


0.8 


Many absorbed AGNs 


P(SF)>50% 


25 


0.60 ± 0.19 


2.20 ± 0.54 


6.3 


1.49 ± 0.81 


2.7 


-0.19 


1.2 


Star-forming galaxies 


















and some absorbed AGNs 


P(AGN)>30% 


34 


0.72 ± 0.17 


1.31 ± 0.45 


1.1 


1.67 ± 0.69 


3.6 


0.12 


0.8 


Many absorbed AGNs 


P(SF)>70% 


13 


0.54 ± 0.15 


2.78 ± 0.77 


5.8 


1.03 ± 1.10 


1.4 


-0.46 


1.7 


Mostly star-forming galaxies 


30%<P(AGN)<50% 


12 


0.77 ± 0.14 


1.56 ± 0.77 


3.1 


2.1 ± 1.2 


2.6 


0.15 


0.6 


Many absorbed AGNs 


(ii) Absorbed AGN Selection (based 


on T = L 2 _ 


10kev/L 


[O iii] A5007* see 










log(T)<0.25 e 


13 


0.71 ± 0.19 


1.72 ± 0.68 


3.9 


2.7 ± 1.0 


4.0 


0.22 


0.1 


Very absorbed AGNs 



a Number of galaxies in the X-ray stack. 

Count rate in units of count s — 1 . 

c Hardness Ratio = (H-S)/(H+S), where H and S are the X-ray counts in the hard (2-8 kcV) and soft (0.5-2 kcV) bands. 
d Effective Photon Index. 

c Also required P(AGN)>30% and log(M*[M©]) > 10.2 as selection criteria. 



